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ABSTRACT 

Results on the properties of warm molecular hydrogen in 57 normal galaxies 
are derived from measurements of H2 rotational transitions in the mid-infrared, 
obtained as part of the Spitzer Infrared Nearby Galaxies Survey (SINGS). This 



CLc study extends previous extragalactic surveys of emission lines of H2, the most 

O ■ abundant constituent of the molecular interstellar medium, to fainter and more 

^ • common systems (LpiR = 10^ to 6 x 10^'^ L0) of all morphological and nuclear 

, types. In our sensitive integral-field observations covering kiloparsec-scale areas, 

the 17/im S(l) transition is securely detected in the nuclear regions of 86% of 

"p^ • SINGS galaxies with stellar masses above 10^'^ M©. The derived column densi- 

^ . ties of warm H2 (with temperatures above ~ 100 K), even though averaged over 

Cf^ . large areas, are commensurate with values observed in resolved photodissociation 

^ . regions; the median of the sample is 3 x lO^^cm"^. They amount to a signifi- 

l> 
O . 

. and more than 30%. The power emitted in the sum of the three lowest-energy 

transitions is on average 30% of the power emitted in the bright [Sill] cooling 
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line (34.8 /iin), and represents about 4 x 10~^ of the total infrared power within 
the same area for star-forming galaxies, which is consistent with excitation in 
photodissociation regions. The fact that the H2 line intensities scale tightly with 
the emission in the aromatic bands, even though the average radiation field in- 
tensity within the same area varies by a factor ten, can also be understood if 
both tracers originate predominantly in photodissociation regions, either dense 
or diffuse. A large fraction of the 25 targets classified as LINERs or Seyferts, 
however, strongly depart from the rest of the sample, in having warmer H2 in 
the excited states, smaller mass fractions of H2 in the warm phase, and an excess 
of power emitted in the rotational transitions with respect to aromatic bands, 
the total infrared emission and the [Sill] line. We propose a threshold in H2 to 
aromatic band power ratios, allowing the identification of low-luminosity AGNs 
by an excess H2 excitation. A dominant contribution from shock heating is fa- 
vored in these objects. Finally, we detect in nearly half the star-forming targets, 
in particular in low-density central regions of late-type galaxies, non-equilibrium 
ortho to para ratios, consistent with the effects of pumping by far-ultraviolet 
photons combined with incomplete ortho-para thermalization by collisions, or 
possibly non-equilibrium photodissociation fronts advancing into cold gas. 

Subject headings: galaxies: ISM - infrared: galaxies - infrared: ISM - ISM: lines 
and bands - ISM: molecules - surveys 



1. Introduction 



Rotational transitions of molecular hydrogen, lying in the mid-infrared range between 
5 and 30 /xm (Tabled]), provide measurements of the mass and temperature distribution of 
the bulk of the warm molecular gas phase, at temperatures of ~ 100 to 1000 K. By contrast, 
the transitions from higher vibrational levels in the near-infrared, better studied because 
more easily observable from the ground, arise from gas at apparent excitation temperatures 
of more than 1000 K located in a thin laye r of molecular clouds, hosting a negligible fr action 
(of the order of 10~^) of the total H2 mass (IBlack &: Dalgarndll976l : iBurton et al.lll992l ). The 
rotational lines are thus more appropriate tracers of molecular gas being exposed to mod- 
erate heating, and arise from a much larger volume fraction of the mo lecular clouds. They 
const itute one of the most important coolants of warm molecular gas (INeufeld fc Kaufman 
1993h . 



H2 emis sion lines have been detected in a wide array of sources, including outflows from 
young stars f Gautier et al.lE976 : Bally fc Lane 1982 ). photodissociation regions ( Gatley et al. 
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19871 : iTanaka et al. Il989l). planetary nebulae ( 



pernova remnants ( Tref ^ra IQTqI: Burton et al 



(iThompson et al 



1978 



lisions (IHerbst et al 



reffers et al.lll976l : iBeckwith et al.lll978l ). su- 



19891 ) , large regions at the centers of galaxies 



Gatlev et al.lll984) and extranuclear large- scale shocks in galaxy col- 



199G; 



Sugai et al.l 119971 : lAppleton et al.ll2006l ). The possible excitation 



mechanisms are accordingly varied. In normal galaxies, the major excitation source is ex- 
pected to be the far-ultraviolet radiatio n of massive stars in photod issociation regions, with 
photon energies between 6 and 13.6 eV (IHoUenbach fc TielensI 119971 . and references therein). 
H2 molecules can be pumped by FUV photons into electronically-excited states, followed 
by fluorescence and radiative cascade through the vibration-rotation levels of the ground 
electronic state. Pure fluorescent spectra are produced only if the cascade is not signifi- 
cantly altered by collisions with hydrogen atoms and molecules; if the critical densities for 
collisional deexcitation are exceeded, a portion of the pump energy is converted to heat 
by collisions, and the lowest rotational levels are populated by collisions and thermalized. 
Pure fluorescence is thus much more likely in the vibrational transitions, that have high 
critical densities, than in the pure rotational transitions considered here, with critical densi- 
ties below a few lO^cm"'^ for S(0) to S(3). Additionally, FUV photons can be absorbed by 
dust grains, followed by the ejection of photoelectrons that heat the gas. This also results 
in the thermal excitation of the low-energy levels of H2 by collisions with the warm gas. 
Besides the radiation of massive stars, a second important source of excitation is shocks, 
in molecu lar outflows, supernova remnants or cloud c ollisions in a disturbed gravitational 
potential (jShuU fc Hollenbaclj|l978l : iDraine et al.lll983l ). In addition to the above processes. 



X-rays produced in active nuclei or in supernova remnant shocks can partially ionize and 
heat the gas over large column densities, leadi ng to H2 excitation by collisions with hydrogen 



atoms and molecules , and with fast electrons (ILepp fc McCraylll983l : iDraine fc Woodslll991 



Maloney et al.l 1 19961 ). Finally, H2 molecules can be formed directly into excited states. 



Surveys of molecular hydrogen line emission in galaxies have been so far mostly restricted 
to starbursts, active galactic nuclei and ultraluminous systems, and have been performed 
mostly in the near-infrared, targetting vibration-rotation lines that arise from upper levels 
with much higher excitation energies than the mid-infrared lines. It has been speculated 
that the major source of H2 heating in star-forming galactic nuclei was shocks in super- 
nova remnants, based on comparison of the luminosity of some vibrati on-rotation H2 lines 
with a limited number of Galactic templates and with shock models ( iMoorwood fc Oliva 



19881 : iMouri et al.lll990l ). However, scaling individual templates to the integrated emission 



of galaxies has large inherent uncertainties, and the near-infrared line ratios most often 
used to discriminate between thermal and non-thermal emission are not always sufficient 
to distinguish between s hocks and fluorescent excitat io n followed by collisio nal deexcitation 



in high-density regions (ISternberg &: Dalgarnd Il989l ). iPuxley et al.l (119881 ) surveyed star- 
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burst galaxies in several vibration-rotation lines, and found that the dominant excitation 



mechanism was pu mping by the 



ar-ultraviolet radiation of massive stars, rather than col- 



lisional excitation. iDavies et al.l (120031 ) reached the same conclusion for a small sample of 



ultraluminous galaxies, in which the first vibrational level is thermalized by high densities in 
photodissociation regions. Active nuclei (LINERs or Seyferts) can show a n excess of H2 emis- 
sion relative to hydroge n recombination lines and aromatic bands (e.g. iMoorwood fc Oliva 
19881 : iLarkin et al.lll998l ). but the exact nature of the additional source of excitation, namely 
X-ray excitation, fluorescence induced by a non-thermal ultraviolet continu um, or shocks 
induced by dynamical perturbations, is often unclear (e.g. iQuillen et al.lll999l ). It is however 



unlikely that significant H2 e mission could arise f r om interaction be tween molecular clouds 
and jets from Seyfert nuclei (IRotaciuc et al.lll99ll : iKnop et al.ll200ll ). 



The detection of a rotational line of H2 was first reported by lBeck et al.l (119791 ) (the S(2) 
transition at 12.3 /im in Orion) from observations at Las Campanas Observatory. It was soon 
followed by many more ground-based detections, but the majority of data on the rotational 



spectrum of H2 were pr oduced by the SWS instrument on board ISO (e.g. iLutz et al.l 12000 



Rigopoulou et al.ll2002l ). Furthermore, with previous infrared spectroscopic capabilities, ob- 
servations of normal galaxies have proven difficult due to sensitivity limitations, so that our 
current knowledge is mainly extrapolated from studies of very bright objects, maybe not 
representative of the general galaxy population. The purpose of this paper is thus to extend 
previous work to fainter systems than formerly accessible, and to characterize directly the 
generic properties of the warm molecular hydroge n content of nornial ga laxies. The SINGS 
sample (Spitzer Infrared Nearby Galaxies Survey; iKennicutt et al.ll2003l ). covering a broad 
range of infrared luminosities, morphologies and nuclear types, is ideally suited to such a 
pursuit. 

Studies of rotational lines alone, without information on vibrational levels, have very 
limited diagnostic value concerning the source of excitation, because the low critical densities 
of the rotational levels make it likely that they will be thermalized most of the time, and thus 
cannot be used to distinguish between the various heating mechanisms. Because observations 
of vibration-rotation transitions in the near-infrared are still scarce for normal galaxies, and 
because they are typically performed in apertures that are not matched to our observations, 
we did not attempt to include vibrational levels in our analysis. The characterization of 
excitation mechanisms and physical conditions in the gas would greatly benefit from such 
information, but would necessitate an additional dedicated survey. 

The rotational lines are, however, energetically important and can characterize the 
temperature and density conditions of a large mass fraction of the interstellar medium in 
galaxies, i.e. that consisting of warm molecular gas. From a SWS survey of rotational lines 
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in nearby starburst and Seyfert galaxies. IValentijn et al.l (Il996l ) and lRigopoulou et al.l (120021 ) 
obtained mass fraction s of in the warm p hase of several percent. In ultraluminous galaxies 
observed with Spitzer, iHigdon et al.l (120061 ) derive much lower mass fractions of warm gas, 
but the fact that the majority of their sample has only upper limits for the S(0) line makes it 
possible that the temperatures are overestimated (because computed from the S(l) to S(3) 
lines only, whenever S(0) is undetected) and thus the masses of warm H2 underestimated. 

This paper presents observations of warrn molecular hydrogen in nearby galaxies ob- 
tained as part of SINGS (IKennicutt et al.ll2003l ). From these data, we present quantifications 
of the temperatures and column densities of warm H2 encountered in kiloparsec-scale areas, 
mostly nuclear regions, and a comparison of the power emitted in the rotational lines with 
those produced by [Sill] at 34.8 fim, which is the dominant cooling line of normal galaxies in 
the mid-infrared range, and by dust. We emphasize the different properties of star-forming 
regions and nuclei classified as LINERs or Seyferts, and discuss their II2 excitation mecha- 
nisms. The data, analysis methods and observational results are described in Sections [2] to 
[sfl The interpretation of the main findings is presented in Sections [6] and [TJ and a summary 
of the results and conclusions can be found in Section [HI 



2. Data and measurements 
2.1. Targets 



The SINGS sample (IKennicutt et al.ll2003l ). comprising 75 galaxies, is intended to be a 



valuable representative set of local galaxies that are not ultraluminous, and whose moderate 
distances ensure that the properties of the interstellar medium can be studied at relatively 
small spatial scales (a few hundreds of parsecs at the shortest wavelengths). Numerous 
sources with mild starbursts or low-luminosity active nuclei are included. Of this sample, 
we excluded from the present study the objects that were not observed in spectroscopic 
mode because of their very low brightness (DDO 154, Hoi, M81 dwA, MSldwB), or con- 
taining very little dust and nebular emission within the nuclear area mapped by the high 
spectral resolution modules (the quiescent ellipticals NGC584 and NGC 1404, the quiescent 
dwarf galaxies DDO 53, DDO 165, Ho IX and the asymmetric magellanic galaxies NGC 4236, 
NGC 5398, NGC 5408 and IC4710). The dwarf galaxies IC2574 and NGC 5474 were also 
rejected because they lack observations in some of the spectroscopic modules. Of the two 



^For easier comparison to future observations and models, ascii flux tables of all the measured quantities 
are available upon e-mail request. 
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star-forming dwarf galaxies with several extranuclear pointings, Ho II and NGC6822, we re- 
tained only NGC 6822 here; the regions within Ho II are indeed two faint to allow an analysis 
of the H2 excitation diagram, contrary to some regions within NGC 6822. Low-mass galaxies 
with extranuclear pointings will be discussed elsewhere. NGC 3034 (M82) was excluded due 
to the unavailability of nuclear spectroscopy from SINGS, as well as NGC 1377, which con- 
stitutes a galaxy class of i ts own very differen t from the rest of the SINGS sample, and has 



been discussed separately (IRoussel et al.ll2006l ). The sample for H2 measurements comprises 



66 targets in 57 galaxies (Table [2]). The pointings are centered either on the nuclear regions 
(for most targets) or on some bright star-forming complexes (for a few dwarf galaxies and a 
spiral galaxy). Diffuse regions within galactic disks are not covered by the present study. 

The aperture over which we extracted the spectra is the intersection of the various 
areas covered by all four spectroscopic modules. The central position and solid angle of this 
aperture, used to measure all the quantities presented in this paper (line and continuum 
fluxes), is listed for each galaxy in Table [2l In practice, the limiting size is that of the maps 
performed with the high-resolution modules, which were enlarged in a few cases in order to 
cover the emission from a star-forming circumnuclear ring. At the distances of the targets, 
the equivalent linear diameters of the apertures range from 60 pc to 3.8 kpc (distribution 
shown in Fig. [1]), and the median is 900 pc. Although the apertures are in general small 
fractions of the optical extent of the galaxies, the measurements are still averages over very 
large and complex areas. It is expected that a large number of disconnected star formation 
sites, in addition to the nucleus, contribute to the total emission. 



2.2. Broadband imaging 



To estimate flux densities of the dust continuum and of the aromatic bands (also referred 
to as the emission from PAHs, or polycyclic arom atic hydrocarbons ), we used images in the 
3.6/im and 7.9 yum bands of the IRA C camera ([Fazio et al.l 120041 ) . and scan maps in the 
three bands of the MIPS instrument f Rieke et all 2004 ) at effective wavelengths of 24, 71 
and 156 /im. Since in early-type galaxies photospheric emission can make an important 
contribution to 7.9 /xm fluxes, we subtracted an estimate of this component in order to 
obtain pure measurements of aromatic band emission. To this effec t, we scaled 3.6 /xm 
fluxes, assumed to be dominated by stellar emission, as described in iHelou et al.l (120041 ). 
The resulting flux densities are noted -FV.gdust- 



The observing strategy and data reduction are described by iKennicutt et al.l (120031 ). 
The full width at half maximum of the point spread function (PSF) is close to 2" at 7.9 /im, 
6" at 24 fim, 18" at 71 fim and 40" at 156 /im. Flux calibration uncertainties are of the 
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order of 10% in the IRAC bands, and 5%, 10% and 15% in the MIPS 24, 71 and 156 fim 
bands, respectively. To correct for the effects of hght scattering in IRAC arrays, we applied 
to flux densities measured from IRAC maps corrective factors that are appropriate for the 
photometry o f extended sourc es within apertures of arbitrary size (derived by T. Jarrett and 
pubhshed by iDale et al.l 120071 ). For our apertures, the correction factor at 7.9 /im is of the 
order of 10%. 



2.3. Spectroscopic data 



^The targets were observed in mapping mode with the IRS instrument (iHouck et al. 

2004 ). at low spectral resolution between 5 and 38 /im, with the SL and LL slits (A/AA ~ 60 



130) and at high spectral resolution between 10 and 37 /zm, with the SH a nd LH slits (A/AA ~ 
600). The observing strategy is described by lKennicutt et al.l (120031 ) and lSmith et al.l (120041 ). 
The data were pre-processed with the S13 version of the Spitzer Science Center pipeline. 
Pixels with an ab normal responsivity were masked, and spectral cubes were built with the 
Cubism software ([Smith et al.ll2007a] ). The flux calibration was performed as described by 
Roussel et al.l (120061 ). We checked the accuracy of this procedure by systematically comparing 
broadband fluxes from imaging observations and from spectra, and line fluxes from high and 
low spectral resolution spectra, for bright lines that are minimally contaminated by broad 
aromatic features at low resolution (but note that even if the flux calibrations of the different 
modules were in perfect agreement, deviations would be expected from slight misalignment 
between the apertures). We obtain F24(MIPS)/F24(LL) = 1.01 ±0.04 (for targets with F24 > 
0.025 Jy within a diameter of about 50"), F7.9(IRAC)/F7.9(SL) = 0.99 ±0.05 (for targets with 
F7 9 > 0.025 Jy within a diameter of about 30" and accurately determined backgrounds in 
SL maps), F[siii](LL)/F[siii](LH) = 0.96±0.14, F^sm] 34(LL)/F[siii] 34(LH) = 0.92±0.20 and 
-^[siii] i9(LL)/F[siii] i9(SH) = 1.05 ± 0.21 (for targets with line fluxes above 60"). 



2.4. Measurements 

The S(0) to S(3) rotational transitions of H2 (Tabled]) were measured for all targets. In 
addition, we measured the S(4) to S(7) transitions in three galaxies in which these lines are 
bright enough to become detectable at low spectral resolution (see Table [3]). 

In high spectral resolution data, we deflned errors from fluctuations of the pseudo- 
continuum, which was fltted as an affine function of wavelength {F,^ = aA ± b). In SL data, 
errors at each wavelength were estimated from spatial fluctuations of blank flelds within the 
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satellite spectral m aps that are auto matically obtained when the source lies in the other 



half of the slit (see ISmith et al.ll2004l ). Both the fluxes and the errors presented in Table [3] 



were then added linearly for each point of the line profile above the pseudo-continuum. 
The line profiles were constrained to have a width compatible with the spectral resolution, 
since the latter is sufficiently low that no line is resolved. Fig. [2] shows the line spectra 
for the representative galaxies NGC1097, NGC6946, NGC7552, NGC1266, NGC4569 and 
NGC4579. 

The S(l) line is usually the brig htest. Of the non -dwarf galaxies of the SINGS sample 
(with stellar masses estimated as by iLee et al.l (120061 ) above 10^'^ Mq), the nuclear regions 
of 86% are securely detected in the S(l) line, with fluxes above three times the measured 
error. The other 14% are either ellipticals of the LINER type, or late- type spirals (Sc-Sd). 



There are two galaxies in common between this sample and th at of iRigopoulou et al. 

J2002h . namely NGC 7552 and NGC 6946, the latter from the study of lValentiin et al.Nl996h . 
For both, our aperture is larger than the beam of ISO-SWS, which covered an area of 280 
to 380 arcsec^. For the lines that were detected with SWS, we obtain fluxes that are higher 
by factors of 2.3 (S(l) in NGC 6946), 5.6 (S(0) in NGC 6946) and 1.1 (S(l) and S(3) in 
NGC 7552). The exact placement of the ISO-SWS beam is not known. For NGC 6946, 
given this uncertainty, it is conceivable that the H2 emission be t wice as bright in our 8 00 
arcsec^ aperture as in the SWS aperture; but the S(0) line flux of IValentijn et al.l (Il996l ) is 
inconsistent with our data. 

For this study, we estimate total infrared fluxes (TIR) between 3 and 11 00 /xm, deflned 
as a l inear combination of 24, 71 and 156 /im flux densities. The formula of iDale fc Helou 
J2002h is used here , and we have checked that replacing it with the more recent prescription by 
Draine fc Lil (120071 ) does not change the following results in any appreciable way. The infrared 
fluxes are measured within the same area as the other quantities for direct comparison. The 
PSF width at 156 /im is however much larger than the size of our spectroscopic aperture, so 
that some extrapolation is needed. We flrst measure MIPS fluxes within the larger aperture 
used to compare total infrared fluxes with line fluxes measured in the LL module. Then, 
we scale these fluxes by the ratio of F24 measured in the small aperture to F24 measured in 
the larger aperture, which is equivalent to assuming that the spectral energy distribution 
does not change from an area of ~ 300arcsec^ to an area of ~ 2000arcsec^. The associated 
errors are however expected to be small compared with the dynamic range of the quantities 
discussed in Section [51 Simulations of the overestimation of the far-infrared fluxes caused 
by the extrapolation, using a simple model of a point-source starburst (with the spectral 
energy distribution of Mrk 33) superposed on quasi- uniform emission from low radiation fleld 
intensity regions (with the colors of the central regions of NGC 24 or NGC 2403) , indicate that 
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the effect should be in most cases of the order of 20% (when the starburst and quiescent 
components contribute equally at 156 yum), and in extreme cases reach a maximum of a 
factor 2 (when the quiescent component dominates). ISmith et al.l (l2007bl ) reached a similar 
conclusion (see their Section 3.2). 



3. Excitation diagrams 

Excitation diagrams provide a convenient visualization of the distribution of level pop- 
ulations and allow first constraints on the excitation mechanisms (thermal or non-thermal) 
that can produce this distribution. They represent the column density in the upper level of 
each observed transition A^'u, normalized by its statistical weight ^^u, as a function of the up- 
per level energy E^. The flux of a transition can be written as F = hu A f2/(47r), where A 
is the spontaneous emission probability, hu is the transition energy and Q is the beam solid 
angle. In the assumption of local thermodynamic equilibrium, the total column density A^tot 
can be derived from = g^^ A^tot exp(— i?u / (kT)) / Z(T), where = {21 + 1) (2 J + 1) 
is the statistical weight (with the spin number / = for even J or para transitions, and 
/ = 1 for odd J or ortho transitions), and Z(T) ~ 0.0247 T / (1 - exp(-6000 K / T)) is the 



partition function (IHerbst et al.lll996l ). valid for T > 40K 



The apparent excitation temperature can then be derived from each pair of transitions 

by: 

kT = {E^2 - ^ui) / HN^i/N^2 X gu2/9ui) (1) 

with A'"ui/A^u2 = Fi/ F2 X A2/A1 X A1/A2. Since both radiative decay and collisions with H2 
change the rotational number J by an even number, the ortho and para states are largely 
decoupled and should in principle be dealt with independently. 



3.1. Ortho-para thermalization and departures therefrom 



As emphasized by iBurton et al 



(I1992I ). the lower rotational levels of H2 will be in 
coUisional equilibrium over a wide range of conditions, because their critical densities are 
low. Figure [3] shows the critical densities of all the rotational transitions observable with 
the IRS instrument, as a function of temperature, compu ted using the func t ional form for 
the collisional de-excitation rate coefficient by H2 g i ven bylShuU &: BeckwithI (119821 ) and the 
transition probabilities given by lBlack fc Dalgarnd (119761 ). The derived critical densities for 
ea ch line are about an orde r of magnitude lower than those for collisions with H computed 



by iMandy fc Martini (119931 ). the comparison being made at 600 K, since iMandy fc Martin 
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( 119931 ) provide results only for high temperatures. 



The integrated emission from warm H2 in star- forming galaxies is likely to come pre- 
dominantly from the d ensest photodissoci a tion regions (PDRs) w ithin the beam, with den- 



sities above 10^ cm ^ (IBurton et al.l Il992l : iKaufman et al.ll2006l ). in which case the lowest 



rotational levels will be thermalized. Observations of starburst galaxies w ith ISO-SWS 



(IRigopoulou et al.l |2002| ) as well as ultraluminous galaxies with Spitzer-IRS (IHigdon et al. 



20061 ) are indeed consistent with this expectation. At first sight, the same applies to the 



galaxies studied here. 

However, some of the excitation diagrams show departures from thermalization of ortho 
levels with para levels, in the sense that the apparent temperatures derived from each pair of 
transitions of consecutive rotational number are not monotonic as a function of upper level 
energy. Clear examples are NGC 1266 (T(SO - SI) = (201±45) K, T(S1 - S2) = (465±34) K 
andT(S2-S3) = (347±18)K); NGC4254 ((162±9)K, (358±59)Kand (259±38) K); and 
NGC 4631 ((127±8) K, (342±39) K and (268±25) K). Such deviations from thermalization 
can be explained by an ortho to para density ratio in the excited states apparently different 
from the equilibrium value. We have 

QpjD ^ OffihighT Eo (2/0+1) (2J0+I) cy.p{-Eo I (kT)) 
3 Ep(2/p+l) (2Jj,+l) cxpC-Ep / (fcT)) 

_ nPR Eo(2Jo+l) cxp(-i;o / (fcT)) fr,^ 

^^^highT J2j,{2Jp+l) cKpi-Ep I (kT)) 

where the subscripts o and p designate ortho and para levels respectively (Jp = and Jo = 1). 
OPR highT, equal to the actual ortho to para ratio {OPR) in the high-temperature limit, 
expresses deviations from local thermodynamic equilibrium (LTE) if it differs from three. It 
may be called the effective nuclear spin degeneracy ratio, but will hereafter be called the ortho 
to para ratio for con venience. In LTE, OPR ~ 2 for T ~ 100 K and OPR ~ 3 for T > 200 K 



(iBurton et al.lll992l ). but OPi? highT = 3 at all temperatures. Although OPi? highT < 3 may 
be inferred for the excited states (J > 2), this does not imply that the ortho to para ratio 
of the bulk of the gas in the J = 1 and J = states be out of LTE. In the following, LTE 
will refer more particularly to the equilibrium between the ortho and para levels, and not of 
the ortho levels or para levels separately. Extinction effects are discussed in Section 13.31 and 
the interpretation of OPP highT values is postponed to Section [61 

To derive temperatures and column densities, we first determine whether the excitation 
diagram is compatible or not with LTE by inserting explicitly the factor OPP high t/3 in 
the equations for column densities of the ortho levels, and deriving temperatures from each 
pair of consecutive transitions as a function of OPP highT, to verify whether these conditions 
are satisfied: T(SO - SI) < T(SO - S2) < T(S1 - S2) < T(S1 - S3) < T(S2 - S3), since in 
gas with a distribution of temperatures, ratios of transitions with low-energy upper levels 
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always probe lower excitation temperatures than ratios of transitions with higher-energy 
upper levels. T(SO — S2) and T(S1 — S3) are independent of OPi?highT and determined 
directly from the observed fluxes, but T(SO — SI), T(S1 — S2) and T(S2 — S3) depend on 
OPi?highT- For each pair (p, o) = (0,1), (2,1) and (2,3), we have: 

kT{Sp - So) = {E^o- p) / ln(OPi?highT x R) . (3) 

with R = Ep/Fo X Ao/Ap x Xp/Xo x (2 Jo + 1) / (2Jp + 1). Figure H] shows the corresponding 
diagram for two galaxies. In case the above condition on the temperatures is satisfied for 
OPR highT = 3, as illustrated for NGC3198, we fix OPi? highT = 3 ; in the opposite case, 
illustrated by NGC4631, we fit OPP highT as explained below. The excitation diagrams of 
all the galaxies, with fits overlaid, are shown in Fig. [51 

3.2. Temperatures and column densities 

Since in all cases the excitation diagrams indicate that a single temperature does not 
fit all the line fluxes, we assume that the H2 emission is the sum of two discrete components 
of different temperatures, which is enough to reproduce accurately the observed fluxes. In 
the general case of OPRughT = 3, we perform a least-squares fit of the excitation diagram 
to determine the parameters of the two discrete components (the lower temperature Ti, the 
upper temperature T2, and their mass fraction) and the normalization by the total column 
density. The results are listed in Table HI 

When the gas is at a range of temperatures, it is in practice impossible to lift the de- 
generacy between mass and temperature from the lowest-energy levels. Since the column 
density has a very steep dependence on Ti, we adopt two different procedures to fit the exci- 
tation diagrams and ascertain the amplitude of the uncertainties caused by this degeneracy. 
In the first case, we constrain Ti to exceed the value for which the column density is 20% 
higher than the nominal density derived from T(SO — SI). In the second case, we leave Ti 
unconstrained. In the following, both approaches will be retained when discussing results 
that depend on Ti. 

For the results not to be biased by systematic sensitivity differences at the wavelengths 
of the H2 transitions, we also replace the measured errors by a uniform weight. When 
OPPhighT is allowed to be fitted, we fix T2 at 1.3 x T(S1 — S3) in the constrained- Ti fits, 
which was chosen from the median value of T2 in galaxies with OPPhighT = 3. In free-Ti fits 
with OPPhighT = 3, the distribution of T2/T(S1 — S3) is large, with a tail of high values; 
therefore, T2 is first fixed at the median value, 1.14xT(Sl — S3), and then at 1.5xT(Sl — S3), 
to probe the full range of most likely values. Finally, when one flux is an upper limit, we fix 
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both Ti at 0.98 x T(SO — SI) (which increases the total column density by a maximum of 
~ 20% with respect to that obtained with Ti = T(SO — SI) but allows a small contribution 
from hotter gas to the S(0) and S(l) lines), and T2 as above. 

For the three galaxies from which more transitions, up to S(7), could be measured, 
the procedure is the same except that a third component has to be added. The additional 
parameters are T3 and the mass fraction of the second component, and T2 is fixed at 400 K. 

Several galaxies barely satisfy the criterion on temperatures to have OPi?highT = 3, 
with T(S1 - S2) > 0.95 X T(S1 - S3) and T(S2 - S3) < 1.05 x T(S1 - S3). When Ti is 
constrained, the quality of their fits can be improved by allowing OPi? high t to vary. For 
these objects, we provide results with OPi? high t < 3. Allowing OPP high t to be smaller 
than the equilibrium value has the indirect consequence that the derived column densities 
are smaller. The amplitude of this effect is indicated in Table Si Similarly, for NGC 1705 
and NGC 4552, we provide results with OPPhighT < 3 and indicate the change in column 
density with respect to OPPhighT = 3, because although the S(2) transition being an upper 
limit prevents any reliable determination of OPPhighT, the Ti temperatures derived with 
OPPhighT = 3 are the two lowest of the whole sample, raising the suspicion that they might 
be artifacts of the constraint on OPPhighT- We also consider OPPhighT < 3 more likely for 
these galaxies in view of the dependence of OPPhighT on H2 brightness, discussed later in 
Section [61 

The median Ti temperature is 154 K when the fits are constrained (ranging between 
97 and 300 K); when no constraint is applied, the median Ti is 118 K with T2 = 1.14 x 
T(S1 — S3), and 161 K with T2 = 1.5 x r(Sl — S3). The total column densities that we 
obtained, averaged over kiloparsec-scale regions in galactic centers, range between 10^^ and 
2 X 10^^ cm~^ (for constrained-Ti fits), or 2 x lO^^cm"^ (for free-Ti fits), and their medians 
are respectively 3 x 10^°cm~^ and 5-6 x 10^°cm~^ (Fig. E])- This can be compared with 
typical column densities of resolved photodissociation regions in the Milky Way. In the Orion 
Bar, column densitie s of H2 warmer tha r i 400 K, derived fr om rotational lines, lie between 



10^° and 10^^ cm ^ (jParmar et al.l Il99ll : lAUers et al.l l2005l ) . Note that because the Orion 



Bar is observed nearly edge-on, an equivale nt FDR seen face-on would have lower column 
densities. In NGC 7023, iFuente et al.l (119991 ) derived a total column density of 5 x 10^° cm~^ 



for H2 warmer than 300 K. Thus, if the H2 emission in our targets comes from similar 
photodissociation regions, they must occupy in general a very large fraction of the observing 
beam, assuming that they do not overlap on the line of sight. 

Figure E] also shows a clear dependence of the local (nuclear) column density of warm H2 
on the total stellar mass of the host galaxy. The stellar mass and the infrared luminosity being 
correlated for star-forming galaxies, there is a similar dependence on far-infrared luminosities. 
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To first order, the column density of warm H2 shows the same behavior as tracers of molecular 
gas and star formation rate densities, which suggests that the primary source of H2 heating 
is the star formation activity in non-AGN galaxies, and the nuclear regions respond to the 
global mass and luminosity. LINER and Sy nuclei do not follow the correlation shown by 
star-forming regions, and tend to have smaller column densities of warm H2. The differences 
in terms of energy output and excitation mechanisms will be studied in more detail in 
Sections [5l [6] and [71 Since the few extranuclear regions and dwarf galaxies included in the 
sample do not distinguish themselves from the other star-forming targets in any obvious way, 
here and in the following, they are not discussed as separate categories. 



3.3. Optical depth toward H2 



Consistent with the negligible optical dep ths inferred from th e silicate absorption bands 
at 10 /zm and 18 /im in most SINGS galaxies (ISmith et al.ll2007bl'). that support the modest 



values of nebular extinction derived from the Balmer decrement (IDale et al. 



20061 ) , we assume 



zero extinction both in the lines and in the dust con tinuum for al l the tar gets. In eight galactic 
centers among the SINGS sample (included here), ISmith et al.l (l2007bl ) obtained a better fit 
in their decomposition of the low spectral resolution spectra by including a finite optical 
depth in the silicate bands. We expect the warm H2 component to suffer less extinction, 
on average, than the warm dust continuum, because the two emission sources will not be 
cospatial in general, and the regions of high optical depth will be confined to compact regions, 
probably more concent rated than the regio ns participatin g in H-^ emission (see the striking 
example of NGC 1377; iRoussel et al.ll2006l ). In particular, iHigdon et al.l (120061 ) did not see 
any evidence for significant extinction in the rotational H2 lines of ultraluminous galaxies, 
although these objects are expected to have much higher optical depths than the present 
sample. In the absence of any quantitative constraint on the differential extinction between 
the dust and H2, we do not attempt to correct H2 fluxes for extinction. 



Using the extinction law of iMoneti et al.l (|200ll ). valid for the Galactic center, we 
have A(9.7/im)/Av = 0.15, A{28.2 fim)/A{9.7 fim) = 0.25 and A{17.0 i2m)/A{9.7 j^m) = 
A{12.3 fim)/A{9.7 fim) = 0.46. Even assuming the same optical depth toward the warm 
molecular hydrogen as toward the hot dust, the extinction correction would not change sig- 
nificantly the derived column densities. The extinction is modest at 10 /zm, and therefore 
negligible at 28 /im, the wavelength of the S(0) line which dominates the total column density 
determination. Extinction effects would however depress the S(l) and S(3) line fluxes with 
respect to S(0) and S(2), and could thus artiflcially lower the derived OP-RhighT- In the 
following, we put lower hmits to OPRughT values, when less than 3, derived for the eight 



-14- 



galaxies with non-zero optical depth at 10 /im. 



NGC31 98 is the sample gal axy with the highest optical depth in the silicate feature 
according to ISmith et al.l (l2007bl ). but its excitation diagram shows no sign of attenuation 
of the S(l) and S(3) lines relative to the others, and is consistent with OPRughT = 3 
(Fig. H]). The second most obscured galaxy of the present sa mple is NGC 1266 (it also has 
the highest nebular extinction according to iDale et al.l (120061 ). Ay = 4.1 mag), for which we 
derive OPR^ighT < 3. If this were due to optical depth effects, then the S(3) line at 9.7/im 
should be more attenuated than the S(l) line at 17.0 /im. Since this would be consistent 
with the excitation diagram, we cannot exclude that the apparently low OPRughT value 
be an extinction artifact in at least this galaxy. The dissimilar behavior of the two galaxies 
in terms of differential extinction between H2 and the dust could then arise from different 
excitation mechanisms and geometries: whereas in the nuclear regions of NGC 3198, classified 
as purely HII, the H2 emission is presumably distributed over a large volume, the H2 emission 
in the LINER nucleus of NGC 1266 may be much more compact, and not produced by star 
formation processes (see Section [5]). For 13 galaxies with ne gligible silicate extinction in 
the spectral decomposition performed by ISmith et al.l (l2007bl ). the excitation diagrams do 
imply OP-RhighT < 3, whether a constraint on the lower temperature Ti is applied or not. 
In addition, of the 6 galaxies found to have non-zero silicate extinction and OPi?highT < 3, 
three would require t{B.2) > in order to obtain OPR^ighT = 3 after extinction correction 
(by > 25% for NGC 1266, by a factor > 6 for NGC 4631 and by a factor > 3.5 for NGC 5866). 
The three others (NGC 1482, 4536 and 6946) would require either t{B.2) > Tsu, or very low 
Ti temperatures (< 100 K). Since it is unlikely that the optical depth toward II2 be higher 
than toward the dust continuum, we conclude that our finding, discussed in Section |6l is 
robust against extinction effects. 



4. Mass fraction in the warm phase 

In order to estimate the fraction of molecular hydrogen that is heated to temperatures 
above ~ 100 K, we searched the literature for observed intensities of the 2.6mm CO(l-O) 
line within a beam comparable to the solid angle of our observations. Table [5] summa- 
rizes the adopted data. The column density of cold H2 as given here is derived from CO 
velocity-integrated intensities on the main-beam temperature scale, assuming a uniform con- 
version factor of C O(l-O) intensities to H2 column densities of 2.3 x 10^° cm~^/(Kkms~^) 



(jStrong et al.lll988l ). We derived aperture corrections to the CO intensities by projecting on 
a map both the IRS beam and the CO beam. We did not use any deconvolution technique. 
Whenever possible, a map from the BIMA SONG interferometric survey, including the zero- 
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spacing total intensity (iHelfer et al.ll2003l ). was used. Otherwise, we used instead the 7.9 /xm 
map and assumed the spatial distributions of aromatic bands in emission and CO (1-0) line 
emission to be similar at the large spatial scales corresponding to our apertures. This can 
be justified qualitatively by the association of dust with molecular gas and the Schmidt law 
(for a recent study of the spatially-resolved Schmidt law, see Kennicutt et al. 2007, in prepa- 
ration). The applied correction factors are listed in Table [51 In some cases, there are several 
available measurements all giving consistent estimates to within 30% ; the corresponding 
unused references are given within parentheses. 

There are two major sources of uncertainty in this comparison. The first one is inherent 
to the difficulty of matching the physical area covered by the IRS integral-field measure- 
ments, from single-dish or aperture-synthesis measurements within a different beam. The 
second dominant source of uncertainty comes from t he conversion factor of CO intensities 



to II2 masses, assumed uniform here. The result of IStrong et al.l (119881) is deri v ed fro m a 



comparison of Galactic 7-ray emission with CO and HI emission. iDame et al.l (120011 ) ob 



tained a consistent conversion factor by extrapolating the gas-to-dust mass ratio measured 
from HI and far-infrared emission, in areas devoid of CO emission, to molecular clouds. Both 
methods provide an estimate of the total H2 column density, including the warm gas as well 
as the cold gas, for molecular clouds under similar average physical conditions as Galac- 
tic clouds. Note however that conversion factors both significantly lower and significantly 



higher have been derived for normal galaxies. For instance, the recent study of iDraine et al. 
( 120071 ) favors an average value of 4 x lO^'' cm~^/ (Kkms^^), based on global gas-to-dust mass 



ratios in the SINGS sample. In addition, the ratio of H2 column density to CO intensity 
can vary by at least a factor two, d epending on the physical conditions of the regions emit- 



ting in CO (iMaloney fc Blackl 119881 ) . even though our observing aperture is large enough to 
cover a large number of molecular clouds and dilute some of the dispersion in their physical 
properties. In particular, the conversion factor is expected to be lower for compact and 
actively star-forming regions than for more diffuse and more quiescent regions. We discount 
here variations due to metal abundance, since we could find CO measurements for only two 
low-metaUicity targets (NGC2915 and NGC6822_A). 

Figure [7] shows the mass fraction of molecular hydrogen in the warm phase (T > Ti ^ 
100 K) as a function of the minimum temperature of the warm component, as determined by 
the lowest-energy rotational H2 lines. The nuclei classified as star-forming have a relatively 
narrow range of lower temperatures (Ti = 144±24K for 31 nuclei, with or without CO data, 
from the constrained fits). However, for nuclear regions classified as LINERs or Seyferts, the 
spread in temperatures is higher (Ti = 180 ±45K for 25 nuclei). No statistically-significant 
difference exists between the 18 LINER and 7 Sy nuclei. 
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A clear anticorrelation exists between the two quantities plotted (partly the result of the 
degeneracy between temperatures and column densities), which remains intact when restrict- 
ing the sample to those galaxies for which we could find well-matched CO data (i.e. with 
correction factors close to unity and with several consistent measurements). The dynamic 
range in the warm gas mass fraction is much higher than accounted for by the uncertainty on 
the total H2 mass. The uncertainty on the warm H2 mass for individual objects is however 
extremely large, owing to the degeneracy between Ti, often ill-constrained by the data, and 
the column density. The example of NGC4579 is the most striking (see Table Hj). Since its 
rotational levels up to J=5 are close to thermal equilibrium (at a single temperature of the 
order of 300-400 K), such a component at 70 K as found in the free-Ti fit is unlikely to be real. 
Because the fits where Ti is unconstrained allow mass fractions in the warm phase that are 
sometimes unphysical (for example for NGC2976 and NGC4826), we favor the constrained 
fits as more plausible, but emphasize that the mass distribution at low temperatures is in 
general unconstrained. 

In the case of constrained-Ti fits, it appears that for a small set of nuclear regions 
classified as LINERs or Seyferts, the warm H2 phase consists only of a very small fraction 
of the total mass, but heated to higher temperatures than in regions classified as purely 
star-forming. This behavior arises naturally if normal photodissociation region excitation is 
missing, and if the hotter gas is located in a thin layer of molecular clouds, or has a small 
filling factor. In the case of free-Ti fits, only NGC1316 (Fornax A) remains robustly in the 
part of the diagram with high Ti and mass fraction below 3%. The average temperature 
is however still higher for LINERs and Seyferts than for HII nuclei, and the average mass 
fraction in the warm phase likewise lower. The reason for this difference will be further 
discussed in Section [TJ addressing the excitation mechanisms. 

5. Comparison of the powers emitted by warm H2, [Sill] and dust in 

star-forming regions 

In order to empirically quantify the importance of the H2 rotational lines in cooling 
the interstellar medium of normal galaxies, and to put constraints on the possible excitation 
mechanisms of H2, discussed in more detail in Section [6l we examine power ratios of H2 to 
other tracers of the warm interstellar medium extracted from the same observations. The 
results presented here are independent of any fits to the excitation diagrams. Only the HII 
nuclei and complexes are considered, LINER and Sy nuclei being separately discussed in 
Section [71 Since the bulk of warm H2, at the lowest rotational temperatures, emits mostly in 
the S(0) to S(2) lines, whereas the S(3) line emission has a noticeably higher contribution from 
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hotter H2, probably indicating more mixed excitation sources (anticipating the discussion of 
excitation mechanisms, see Section [6TT]) . we choose, as the most useful quantification of H2 
rotational emission in star-forming targets, the sum of the S(0) to S(2) lines. 



5.1. Total infrared emission 



In photodissociation regions, almost all the far-ultraviolet power from massive stars 
that does not escape is absorbed by dust and converted to infrared continuum radiation, 
or is absorbed by H2. Only a very small fraction of the power absorbed by dust, of the 
order of 1%, is converted to photoelectrons that heat the gas, and emerges as infrared lines 
( iTielens fc Hollenbaclj|l985l ). The dominant gas coolants are the [Ol] and [Cll] lines at 63 fim 
and 158 /xm, but mid-infrared lines, in particular [Sill] at 34.8 /xm and the H2 rotational lines, 
are also energetically significant. Although the transition rate coefficients of H2 are low and 
the excitation energies relatively high, H2 molecules are dominant in number. 

The observed ratios of the power emitted in the sum of the S(0) to S(2) lines to the 
total dust power emitted in the infrared (TIR; see Section [2^ range between 2.5 x 10~^ and 
7.5 X 10~^ for nuclear regions that are not classified as LINERs or Seyferts ( Fig. [8^) . These 



ratios are in agreement with predictions of the photodissociation models of iKaufman et al. 



(I2OO6I ) for a wide variety of radiation field intensities Gq and hydrogen densities n, but a 
relatively narrow range of Gq/u ratios, approximately between 0.1 and 1 with Go in units 
of 1.6 X lO^'^ergs"^ cm~^ and n in units of cm~^. Note that models predict the ratio of 
the H2 line power to the far-ultraviolet (FUV) power (for photon energies between 6 and 
13.6 eV), rather than the total infrared power. Since the intrinsic FUV fiux heating the 
photodissociation regions is unknown, the comparison between observations and models is 
here made by assuming an exact conversion of FUV photons to infrared photons. The 
fraction of dust heating provided by non-FUV photons can however be significant. Allowing 
for this effect would redu ce the derived Gq/ti ratios. The H2 rotational line fiuxes predicted 
bv lKaufman et al.l (12006) are near Iv an order of magnitude higher than those from the older 



models of Burton et al. 



(119921 ). because of the inclusion of photoelectric heating by PAHs, a 



better H2 model, and a finer numerical grid near the region of H2 emission. 



The inferred Go/n ratios are lower than the results of iMalhotra et al.l (120011 ). who de- 
rived the physical conditions of an ensemble of bright star-forming galaxies from the [Cll] 
and [Ol] lines. They found Go/n ratios between about 0.5 and 6, i.e. on average 5 times 
higher than those indicated here by the rotational H2 lines. A possible explanation is that 
H2 emission comes from cooler and denser regions tha n [ClI] and [Ol], because H 2 exists at 
higher optical depths inside the clouds than C"*" and O (IHoUenbach fc Tielenslll997l ). The dif- 
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ference in physical conditions could thus merely reflect a different spatial origin. Besides the 
different locations within PDRs, th e two s tudies also deal with different regions within galax- 
ies: the targets of iMalhotra et al.l (l200ll ) were selected to have most of their line emission 
encompassed by the ISO-LWS beam of 70", whereas our apertures usually c over small frac- 



tions of the line and dust emitting areas. Alternatively, the observations of iMalhotra et al. 



(120011 ) could reflect intrinsically different physical conditions because their sample contains 
galaxies on average brighter and more active than the sample used here. Their far-infrared 
luminosities (in the definition of iHelou et al.lll988l ) range from 6 x 10'' to 8 x 10^^ Lq, with 
a median of 1.5 x 10^° Lq, whereas the far- infrared luminosities of the present sample range 
from 10'' to 6 X 1 0^'^ Lf^, with a median o f 3 x 10^ Lq. The median Fgo/-Fioo ratio is also higher 
in the sample of IMalhotra et al.l (l200ll ) (0.57) than in our sample ( 0.41), indicating higher 
radiation field intensities on average. The Gq/ti ratios derived by IMalhotra et al.l (1200 ll ) 
however do not display any clear correlation with either infrared luminosity or color. Only 
NGC 1482 and NGC 5713, included in both samples, allow a direct comparison of model 
results (we discard the LINER NGC 1266 because most of its H2 emission is not produced 
by PDRs, as shown in Sect. [7]). For both sources, the H2 line fiuxes indicate consistently 
Go ~ 400 and n ~ 1-2 x 10^ For NGC 1482, Gq is in agreement with one of the two 
models of IMalhotra et al.l (I2OOII), but n is at l east f our times higher. For NGC 5713, Go is 
two times higher that that of IMalhotra et al.l (I2OOII ). and n is at least six times higher. In 
conclusion, we favor differences in spatial origin (both within PDRs and within galaxies) as 
a likely cause for the different model results. 



5.2. [Sill] line emission 

Figure [Ht> shows the ratio of powers emitted in the H2 rotational lines and in the [Sill] 
line. The dispersion in the ratio is very similar to that seen in Fig. [8^, and the [Sill] line 
alone emits more power than the sum of the S(0) to S(3) transitions in HII nuclei. The 
[Sill] line has indeed been found to be the brightest mid-infrared cooling line and to scale 
tightly with the total infrared power both in nuclear and extranuclear regions within the 
SINGS sample galaxies (Helou et al., in preparation), with only a very slight dependence 
on the radiation field intensity. We have on average F(SO — S2)/F([SiII]) = 0.3 (ranging 
between 0.15 and 0.5 for nuclei), and F( [Sill] )/T/i? = 2 x 10"^ Using the [Sill] line as a 
substitute for the total dust emission is advantageous because it is observed at about the 
same angular resolution as the H2 lines, whereas estimating the total infrared power within 
these apertures requires a large extrapolation (because of the large width of the point spread 



function at 70 and 160 /xm), making the uncertainty on H 9/TIR relatively 



power predicted by the photodissociation region model of iKaufman et al.l (120061 ). with the 



ligh The [Sill] 
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same physical conditions as above, is however smaller than observed by a factor greater 
than 3, which implies either that the majority of [Sill] emission comes from HII regions in 
high-met allicity nuclear regions, or that the fraction of silicon incorporated in dust grains is 
smaller than 90%. 

Only the regions B and C in NGC6822 have significantly less [Sill] emission, with 
respect to H2 emission, than the nuclear regions of spiral galaxies. Their H2 emission is 
also slightly overluminous with respect to the aromatic bands (Fig. [Hb). This may not be 
entirely attributable to a metallicity effect, decreasing the abundances of PAHs and silicon, 
since region A (Hubble V) has n ormal flux ratios, and oxygen abundances are quite uniform 
in NGC6822 (jPagel et al.lll980l ). An alternative explanation is that additional excitation of 
H2 may be provided in regions B and C, with respect to region A, by shocks in supernova 
remnants (see the more general discussion in Sect. 17. 2p . To our knowledge, no ind ependant 
evidence exists to test the existence of shocks in these regions. IChandar et al.l (I2OOOI ) obtained 
a normal HII optical line spectrum at the center of NGC 6822_C, but since their beam of 
2.5 arcsec^ is only about 1% of ours, we cannot rule out shock excitation. Finally, given the 
small distance of NGC 6822, the regions covered by the IRS aperture are less than 100 pc 
in size. Greater fluctuations around the average properties are thus not unexpected. At 
present, we are unable to decide which scenario is the most likely. 



5.3. Aromatic bands 



Figure [Hb shows a remarkable constancy of the power ratio of the H2 rotational lines to 
the aromatic bands. Among the measured dust and gas observables, PAH emission provides 
the tightest correlation with H2. Observations of photodissociation regions have shown that 
the emission from aromatic band carriers and from fiuorescently-excited H2 just outside pho- 
toionized regions are nearly cospatial, with H2 sometimes seen to extend slightly deeper into 
mole c ular clouds (ISellgren et al.lll990l : iGraham et al.lll993l : iTielens et al.lll993l : iBrooks et al. 



2000; Habart et al. 



2OO3I ). Cospatiality might be expected since both species can be excited 



by FUV photons. Aromatic band carriers can also be excited by lower-e nergy photons in th e 
ultraviolet and optical, but with smaller absorption cross-sections (see iLi fc Draindl200ll ). 
so that FUV photons will dominate the excitation whenever massive stars are present. H2 is 
however dissociated by FUV photons between 11.3 and 13.6 eV where it is not self-shielded, 
whereas PAHs survive the absorption of these photons. Therefore, in the case of relatively 
dense PDRs (associated with molecular clouds), where coUisional heating is expected to be 
the major origin of the H2 rotational lines, H2 emission should peak at slightly higher optical 
depth than aromatic bands, in the transition layer between atomic and molecular hydrogen. 
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with Ay > 1. In addition, PAHs probably cannot be excited as deep into molecular clouds as 
H2, because at sufficiently high de nsities they will be coagulated onto grain mantles on short 
timescales ( iBoulanger et al.lll990l ). If photodissociation regions dominate the excitation of 
H2, as consistent with the above results, a tight relation between aromatic band emission 
and rotational H2 emission can arise only if the physical conditions in PDRs, especially the 
Go/n ratio, are relatively uniform, because H2 fluxes and PAH fluxes depend in very differ- 
ent ways on these two parameters. The condition of relatively constant Gq/u ratios seems 
verifled in the present sample at least for the average emission wit hin kiloparsec-scale re- 
gions (see above). Based on the modelling of [ClI] and [Ol] emission, iMalhotra et al.l (I2OOII ) 
proposed that a regulation of Gq/u might be achieved at the scale of individual PDRs by 
expanding HII regions in pressure equilibrium with their surrounding PDRs. 



A correlation was previously claimed by iMouri et al.l (Il990l ) between the 3.3 fim band 
and the v=l-0 S(l) line at 2.12 fim for a small sample of star b urst and Seyfert galaxies. T he 
dominant source that they propose for H2 excitation, following iMoorwood fc Olival (119881 ). is 
however not photodissociation region s, but shocks in supernova remnants. Using the shock 
models of Kaufman fc Neufeldl (1l996l ) to estimate the sum of the S(0) to S(2) transitions (up 
to 6% of the mechanical power, assumin g that its totalit y is dis sipated in molecular clouds), 
and the population synthesis model of iLeitherer et al.l (119991 ) to estimate both the total 
mechanical power and the FUV luminosity from continuous star formation with a Salpeter 
initial mass function, shocks alone are in principle able to produce a significant fraction 
of the observed H2 emission, but only if the efficiency of conversion of mechanical power 
into II2 emission is unrealistically high. The rotational line ratios are also inconsistent with 
shock models, which predict higher temperatures (T > 1000 K) except for very low shock 
velocities (in which case the power fraction radiated away by rotational H2 lines is lower). 
If the collective rotational line emission from shocks in supe rnova remnants is similar to 
that observed in individual objects such as 3C391 and IC443 (IReach et al.ll2002l ). then this 
mechanism can provide only a modest fraction of the total H2 emission. In addition, if H2 
emission came predominantly from supernova remnants whereas aromatic bands arise mostly 
in photodissociation regions, the partial deconnection between the two, both temporal and 
spatial, would manifest itself by a large scatter in the observed relation between H2 and PAH 
fluxes for galaxies with diverse star formation histories, which is not observed. 



More recen t lv. iRigopoulou et al.l (120021 ) proposed a relation similar to that presented 
by iMouri et al.l (119901 ). between the 7.7/im aromatic band and the rotational S(l) line in 
starburst galaxies. Figure ^ not only confirms this result for lower-luminosity galaxies, but 
also shows that the dispersion for the whole sample of star-forming nuclei is very small, 
and much smaller with the aromatic bands than with the 24 /xm emission (Fig. [Hli), which 
is dominated by the continuum from transiently-heated very small grains, as well as from 
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big grains in intense radiation fields. The quantification of the average H2 to dust power 
ratios and their dispersions is given in the caption of Fig. [HI The energy coupling between 
aromatic band carriers and H2 strongly suggests that both are excited predominantly in 
photodissociation regions, although they may not come from the exact same layers (at the 
same optical depths within the clouds). We present further analysis in the next section. 

A similar correlation, with a similarly small disp ersion, was observed between the [Cll] 
line and aromatic band emission (IHelou et al.ll200ll ). This relation suggests that aromatic 
band carriers are the source of a major part of gas heating in photodissociation regions, via 
the photoelectric effect, at least at modest radiation field intensities, since [ClI] emission is 
the dominant cooling channel in this case. In the narrow range of physical conditions that 
seem to apply if the emission from HII nuclei is interpreted in the framework of photodis- 
sociation region models (a dynamic range in Gq/u of only a factor 10), then the same link 
between aromatic band carriers and H2 would follow if H2 were heated in relatively dense 
photodissociation regions by the PAHs. Our results, however, suggest that in nearly half 
the star-forming targets, the dominant excitation mechanism of the rotational levels may 
be fluorescence in low-density regions, so that ortho-para thermalization is not achieved by 
collisions (see Sect. E]). If the lines are fluorescently excited, the cause underlying the tight 
relation between H2 and aromatic band emission may be that both are proportional to the 
incident far-ultraviolet flux which excited them. 



5.4. Cirrus clouds versus PDRs with high radiation field intensities 



The tight association between H2 emission and aromatic bands (Fig. [Ht) may be sur- 
prising if one assumes that a significant fraction of aromatic band emission arises from 
diffuse, mostly atomic regions with l ow radiation field intensities. The infr ared emission 
of such clouds is often termed cirrus (iLow et al.l Il984j : iTerebey fc FichI Il986l ) . If this were 
the case, then the scaling of PAH flux with H2 flux could be explained only if a constant 
fraction of the total FUV flux escaped PDRs and were absorbed in the more diffuse inter- 
stellar me dium. We stress that we ad opt here the definition of photodissociation regions 
stated by iHoUenbach fc Tielend (119971 ): these are not restricted to the interfaces between 
bright HII regions and dense molecular clouds, but apply more generally to all the neutral 
interstellar medium illuminated by FUV photons (with energies between 6 and 13.6 eV). 
Figure [H] demonstrates the great difficulty of making the idea of an important contribution 
from the cirrus medium consistent with the data. It shows the flux ratios of H2 to PAH 
emission on the one hand, and 24 /im to PAH emission on the other hand, as a function of 
-P24 ~ ^^24-^24/(^^71 -^71 + ^^156-^156), estimated within the spectroscopic apertures. The quan- 
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tity P24 is closely related to fu>ioo = .f(L dust ',Ur».d > 100)^ d erived from the modelling of 



the global spectral energy distributions by iDraine et al.l (120071 ). which is the fraction of the 



total dust luminosity emitted by regions with radiation field intensities f/rad higher than 100 
times the local average value. The dust luminosity fraction of cirrus clouds, ~ fu<io, can be 
evaluated as 1 — c fu>wo, if one assumes that fu>wo and /i7>io are in constant proportion 
to each other. Fig. [9] shows that while the F24/PAH ratio rises by about one order of mag- 
nitude, the H2/PAH ratio is invariant as a function of P24 or fu>ioo- The results obtained 
by replacing P24 with F71/F156, which has the same physical significance as Ff^o /F^no, a mor e 



traditional indicator of the relative importance of HII regions and cirrus clouds (jHeloulll986l ). 
are identical. 

We conclude that aromatic bands are mostly associated with photodissociation regions 
(illuminated by FUV photons able to provide H2 excitation). In addition, since PAHs are 
excited not only by FUV photons but also by low-energy photons, the observed constancy 
of the H2 to PAH ratio imposes some restrictions on possible variations of the radiation 
field hardness. Assuming that cirrus clouds, i.e. PDRs with low radiation field intensities, 
receive appreciably softer radiation than PDRs with high radiation field intensities, it would 
be difficult to understand how both types of regions could produce similar H2/PAH ratios. 
As a corollary, the hypothesis that cirrus clouds could make a large contribution to H2 and 
PAH emission in our targets, although not definitely ruled out, is not favored. Note that 
the situation may be different in more quiescent parts of galaxies, not probed by the present 
sample, and deserves further investigation. Measurement of H2 line fiuxes in quiescent regions 
is however challenging, because they depend steeply on the Gq and n parameters. 

The above does not preclude a large portion of the H2 and PAH emission to originate 
in relatively diffuse molecular gas. Estimates of the optical depth of the ^^CO(l-O) line, over 
large areas of the G alaxy, indicate th at the total molecular medium comprises a substantial 



diffuse component (jPolk et al.lll988l ). We will see in the next section that our data, for a 
portion of the targets, do support an important contribution from low-density PDRs to the 
total warm H2 emission. 



6. Excitation mechanisms in star-forming regions 

In sources with purely stellar activity, H2 emission is expected to arise in varying pro- 
portions from two main energy sources: the far-ultraviolet radiation of OB stars illuminating 
PDRs; and shocks in supernova remnants or other sources, providing collisional heating. In 
the first case, the excitation can be thermal, by collisions with gas heated by photoelectrons, 
or by inelastic collisions with H2 molecules pumped by FUV photons. The excitation mech- 
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anism can also be fluorescence, followed by radiative cascade to the ground vibrational state. 
Heating by supernova remnants is unlikely to be dominant for two reasons, as we have seen 
in Sect. 15.31 First, at the low observed temperatures dominating the warm H2 mass, heating 
would have to be provided by slow shocks, which are not efficient enough to compete with 
PDR excitation. Second, variations in star formation histories within the sample, which are 
shown by Moustakas et al. (2007, in prep.) to be very large (from population synthesis 
fitting to optical spectra), would produce more scatter than observed in the H2/PAH ratio. 
We conclude that H2 is heated predominantly by PDRs, and this interpretation is supported 
both by energetics arguments and by the close association with aromatic band emission 
(Sect. [5]). In this section, we focus on additional constraints on the physical conditions and 
excitation mechanisms in PDRs (thermal or fluorescent) from the line ratios and excitation 
diagrams. 



6.1. Constraints from the temperature distribution 



We have seen that the H2 to far-infrared ratios are consistent with values of Go/n, 
the ratio of the average radiation field intensity to the hydrogen density in PDRs, between 
about 0.1 and 1. In principle, separate constraints on Go a nd on n can be obtain ed from 
the temperature distribution reflected in the H2 line ratios (IKaufman et al.ll2006l ). Given 
the complexity of the surveyed regions, however, the interpretation of the line ratios by 
comparison with models of a single PDR is severely limited. First, the emission from many 
distinct PDRs, presumably showing a large range of physical conditions, is averaged within 
the beam. Second, even though the total emission in the sum of the S(0) to S(3) lines is 
probably dominated by PDRs, shocks also related to the star formation activity must be 
present, dri ven in particular by supernova remnants, protostellar outflows, or turbulence 
dissipation (IFalgarone et al.ll2005l ). Since these shocks are characterized by higher rotational 
temp eratures than PDRs, they would contribute mostly to the S(3) line (IKaufman fc Neufeld 
19961 ). and in view of the observed line ratios, negligibly to the lower- lying transitions. 



In both cases, the superposition of PDRs of different conditions or of shocks induces 
a spread in temperatures, and as a consequence these are not reproduced by single PDR 
models. We com pared the line ratios o f the star- forming targets with OPi? high t = 3 to 
the predictions of IKaufman et al.l (120061 ). Although the ranges in Gq and n derived from 
S(2)/S(0) have a broad overlap with those derived from S(1)/S(0), the S(3)/S(l) ratios are 
inconsistent, indicating more intense radiation fields and higher densities. The least biased 
tracer of physical conditions in the bulk of the PDRs is thus the S(1)/S(0) ratio (for galaxies 
with OPRughT = 3). It suggests that the average Go varies between about 100 and 5000 
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[in units of 1.6 x 10 ^ergs ^ cm ^) 



and the average n between about 500 and 10^ cm 



6.2. Fluorescent excitation 



It is surprising that in several objects, the S(0) to S(3) transitions are not thermalized, 
as indicated by deviations from an apparent OPi?highT of 3. OPR values that are different 
from the equilibrium value at the temperature of the H2 gas arise naturally from fluorescent 
excitation, if the gas density i s lower than the critical de nsity for ortho-para equilibration 
by collisions with H and H"*" ( ISternberg fc Neufeldl Il999l ). This is because the ultraviolet 
absorption lines have a greater optical depth in the ortho states than in the para states, so 
that the ortho states at a given depth are pumped less. The apparent OPi? high t values only 
apply to the states excited by FUV pumping, and do not imply that the true OPi? high t 
value (for the J=l and J=0 states, where most of the gas resides) be different from 3. 
For a total OPPhighT of 3, the OPPhighT of the excited levels is predicted to be close to 
-\/3 = 1.7. A thorough review of the phenomenon of selective e xcitation and its impl i cation s 
for the interpretation of excitation diagrams was provided by ISternberg fc Neufeldl (119991 ). 
Note that radiative decay in the electronic ground state and most collisional (de)excitations 
always occur at fixed spin number and thus preserve the ortho or para state. Conversion 
from one state to the other can be accomplished by H2 dissociation followed by reformation 
on dust grains, or by reactions with protons and hydrogen atoms in the gas phase. 

Deviations from t hermalization had previously remained unse en for rotational lines of 
extragalactic sources (jRigopoulou et al.l l2002l : iHigdon et al.ll2006l ). In the frame of PDR 
excitation, values OPPhighT < 3 may be interpreted as arising in FUV-pumped gas with 
sufficiently low densities to prevent ortho-para equilibration by collisions. However, because 
this process depends on unknown timescales for gas heating, cooling, dissociation and refor- 
mation, we are unable to quantify in a simple way the implied density conditions. Alterna- 
tively, it is conceivable that the emission comes from initially cold gas, that h as been heated 
by slow shocks and has not had time to reach the equilibrium value of OPR ( iTimmermann 
19981 : IWilgenbus et al.ll2000l ). For HII nuclei, however, we have seen that photodissociation 
region excitation is more likely than shock heating (Sect. [5]). 

Another possibility is that in a fraction of the PDRs within the beam, the photodis- 
sociation front is advancing into cold gas (in LTE at time ti with OPR{ti) ^ 3 and 
OPP high t(^i) = 3), and the recently-heated gas has not yet had enough time to reach 
LTE at time ^2 (with 0PR{t2) ~ O PR(tl) and OPRhs^ hTih) < 3 according to E qu. [2]). 
This is the interpretation favored by iFuente et al.l (119991 ) and iHabart et al.l (120031 ) to ex- 
plain the non-LTE ortho to para ratios observed in the PDRs of NGC 7023 and p Ophiuchi, 
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respectively. In this scenario, the observed portion of the interstellar medium would have 
to contain a much larger fraction of non-equilibrium PDRs in targets with low OPi? high t 
than in targets with OPi?highT = 3, and the underlying reason would be unclear. It is also 
unknown whether the timescales involved in ortho-para equilibration in our sources are long 
enough for this scenario to be viable. 

Figure [TO] shows the derived OPPhighT values (fixed to 3 whenever the temperatures 
derived from each pair of adjacent transitions were compatible with this assumption) as a 
function of the total brightness of the S(0) to S(2) lines. We have here included LINERs 
and Seyferts because they do not display a different behavior in this diagram. Under the 
hypothesis of photodissociation region excitation, and assuming first that sites of star for- 
mation occupy a constant fraction of the observing beam, sources with low H2 brightnesses 
should consist of regions with both low densities and low radiation field intensities, while 
sources with the highest H2 brightnesses should include a greater fraction of high-density, 
high-radiation regions. In this simplified view, low values of OPP high t, indicating that H2 is 
not thermalized by collisions, could be obtained only in the low-brightness sources, as seen 
generally in Fig. [TOl Variations in the beam filling factor by sites of star formation would 
then induce a horizontal scatter, which is indeed very large. Fig. (TU] also shows OPP high t as 
a function of the average surface brightness in the 24 /xm band. The latter quantity incorpo- 
rates a significant contribution from HII regions, and should be dominated by variations in 
radiation field intensity, rather than variations in gas density. The fact that the horizontal 
spread is larger in this diagram than in the diagram involving the H2 brightness supports our 
tentative interpretation in terms of density effects. In addition, OPP high t does not show any 
variation as a function of P71/P156, and only a weak tendency to increase with the quantity 
P24 discussed in the previous section. Both these quantities are indicators of the average 
radiation field intensity; they are observed to be generally correlated with the gas density, 
but only weakly. Although the data do not allow us to truly estimate average densities in 
photodissociation regions, they suggest that in a substantial number of the observed nuclear 
regions, the emission can be dominated by low density gas, relative to well-studied Galactic 
PDRs and starburst galaxies. New modelling is required to quantify the conditions under 
which rotational lines indicate OPP high t < 3 in PDRs. 



The H2 line ratios, compared with the PDR model of iKaufman et al.l (120061 ) . do not 
indicate lower gas densities, on average, in galaxies with OPPhighT < 3 than in galaxies 
with OPPhighT = 3. However, the densities estimated in this way are averages within the 
whole beam, and the densest and warmest regions have a greater weight, because they are 
more luminous in H2; if a large spread in densities exists, with both dense clumps and diffuse 
PDRs, an increase in the diffuse fraction may not be easily detectable in the average density, 
while still leaving an imprint on OPPhighT- We note that the proton density should also 
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play an important role, since ortho-para conversion is effected by collisions with H and H+. 
Whether galaxies with OPi^MghT < 3 have lower proton densities in molecular clouds, due to 
reduced ionization by cosmic rays, is in principle testable with radio continuum observations 
of synchrotron radiation from cosmic ray electrons. The number of sample galaxies with 
adequate data, at sufficiently high angular resolution, is however too small to apply this 
test. 

We have seen in Section 13.31 that extinction effects are unlikely to modify our results 
in a statistical sense. The fitted OPi?highT values are correlated with the brightness of the 
dust emission, in the aromatic bands, the 24 /xm continuum and total infrared emission (not 
shown here but similar to Fig. [TOl) . which is a further indirect argument against extinction 
effects being responsible for low OPi? high t values. If extinction played a significant role, then 
low OPP high T values would be seen preferentially in bright and compact regions, which is 
not the case. 

The galaxies for which we derive the lowest OPP high t values (and among the lowest 
H2 surface brightnesses, regardless of the constraint on the Ti temperature) are NGC337, 
NGC1705 (OPPhighT is a lower limit), NGC2915, NGC4552 (OPPhighT is a lower limit) 
and NGC 7793. The first three are dwarf galaxies, NGC 7793 a very late-type spiral and 
NGC 4552 is a small elliptical galaxy classified as LINER, with the smallest infrared brightness 
of the whole sample. Since smaller gas densities are expected in general in dwarf galaxies 
than in the central regions of massive galaxies, except in blue compact dwarfs such as Mrk33, 
this finding is consistent with our interpretation of small OPPhighT values in terms of low 
density. We may also remark that i n NGC 1705, which is a starburst, no ultracompact 



HII region was detected in the radio (I Johnson et al.ll2003l ): thus no conflict exists with the 



hypothesis that the photodissociation regions in NGC 1705 have low densities. 



7. Excitation mechanisms in LINER and Seyfert nuclei 

A large number of the galaxies classifled as LINERs or Seyferts deviate significantly from 
the relations discussed in Section [5l in having a strong excess of H2 emission with respect to 
all the other tracers used here, not only aromatic bands but also [Sill], the 24 /xm fiux and 
the total infrared emission (Fig. [8]), arguing for an alternative excitation mechanism in these 
galaxies. The average and dispersion of each power ratio are given separately for HII nuclei 
and for LINER and Sy nuclei in the figure caption. Given these results, the quantities by 
which the two categories are most clearly separated are the H2 to aromatic band ratio and 
the H2 to [Sill] ratio. In particular, nuclear regions with P(SO — S2)/P7.9dust > 10~^'^^ are 
likely to be of the LINER or Sy type at the 99% confidence level. 
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We thus define an excess of H2 emission, witli respect to HII nuclei, based on tfie observed 
relation with aromatic bands. We choose as the maximal H2 power associated purely with 
star formation the quantity 10~^'^^ x L7 9dust- Aromatic band carriers are thought to be 



destroyed in intense radiation field enviro nments, such as HII region cores (IGiard et al.lll994f) 



and ionized regions around Seyfert nuclei ( Desert fc Dennefeld 1988 : Voit 19921 ). They would 



however survive where H2 is not dissociated, so that an enhancement of the H2/PAH ratio 
is more likely caused by a genuine excess of H2 emission, in conditions where PAHs are not 
excited, rather than normal H2 emission occurring where PAHs would have been destroyed. 
Furthermore, H2 emission is seen in excess not only with respect to aromatic bands, but 
also with respect to [Sill] and other dust tracers (24 /im and total infrared emission). [Sill] 
emission may be depressed because the ionization state of silicon becomes higher, but the 
dust continuum cannot be suppressed like aromatic bands. Our empirical quantification 
of the H2 excess is intended to extract the part of H2 emission that cannot originate in 
photodissociation regions. The excitation mechanism that is required to account for this 
excess, while not exciting PAHs, is either X-ray irradiation or shock heating. 

We discuss LINER and Sy nuclei as a single category, because no detectable difference 
exists in their H2 properties. This may be expected for several reasons: small-number 
statistics; the fact that the classification of low-luminosity AGNs can be ambiguous as it 
depends on the aperture size in particular; and the fact that the source of H2 excitation 
might not be directly linked to the nuclear activity, as the results discussed below suggest. 



7.1. Heating by X-rays from an active galactic nucleus 

The idea that a nuclear X-ray source may modify the chemistry and excitation of the sur- 
rounding molecular clouds through sufficiently large column densities as to offer a convenient 



diet, in particular, unusua l ratios of tracers of dense molecular gas (ILepp &: Dalgarno 



(LeDD & Dalgarno 


1996; 


8( 


Usero et al. 


2004 


). X- 



ray excitation would also manifest itself in the properties of H2, the most abundant molecule. 
To test the hypothesis that the additional excitation in galaxies showing a significant excess 
of H2 emission with respect to aromatic bands (Figure ^) is predominantly produced by 
nuclear X-rays, we have compiled estimated X-ray fluxes in the 2-10 keV band obtained from 
Chandra observations. The data and references are summarized in Table O The H2 excess 
is shown as a function of the X-ray luminosity of the nucleus in Fig. [TTJ Here, galaxies with 
no H2 excess according to our definition are shown below the dashed line. 



The spread in Figure [TT] is very large. Since the available X-ray measurements do 
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not isolate the hard X-ray component, but include soft emission, a substantial part can 
be thermal emission from supernova remnants, as opposed to power-law emission from the 
activity related to a central supermassive black hole. Many HII nuclei with no H2 excess 
indeed have X-ray luminosities that are comparable to those of LINER and Seyfert nuclei. 
On the other hand, some galaxies, in particular NGC3627, NGC4569 and NGC5866, are 
very luminous in H2 but have only modest X-ray to H2 luminosity ratios compared with other 
H2-excess galaxies. We shall assume a power-law spectrum with a standard photon index 
of —1.8, as adopted by iHo et al.l (120011 ). to extrapolate the total X-ray luminosity to lower 
photon energies. Under this assumption, the luminosity between 0.2 and lOkeV is only two 
times the luminosity between 2 and 10 keV compiled here. Up to 10% of the intrin sic 1-10 keV 
luminosity may emerge in the sum of all H2 transitions (iLepp fc McCraylll983l ). whereas in 
our comparison we have summed only the three rotational lines S(0) to S(2). Therefore, 
pure X-ray excitation of the excess H2 emission should not be energetically possible even in 
objects such as NGC4579 and NGC5195, which have relatively low H2 to X-ray luminosity 
ratios in Fig. [TTl Even if the intrinsic emission in soft X-rays were underestimated by a large 
factor. X-ray excitation would still be unlikely to dominate in most cases. The apparent 
trend of increasing excess H2 luminosities with increasing X-ray luminosities in Fig. [TT] does 
not imply a direct excitation of H2 by X-rays, since it could also be understood in the frame 
of multi-phase shocks, produce d by supernova rern nants or by starburst winds (see below). 
From near-infrared line ratios, iDavies et al.l (120051 ) also reached the conclusion that X-rays 
do not contribute significantly to the excitation of H2 in a small sample of active galactic 
nuclei, within regions smaller than those sampled here (about 100 pc). 

Supernova remnants are another important source of X-rays. More than h alf their me- 
chani cal energy can be converted to X-rays, mostly with energies below ~ 500 eV (jPraine fc Woods 
I991I ). but the X-ray powe r is converted to em ission in the S(0) to S(2) lines with low 
efficiencies, less than 10~^ (jPraine fc Woods! 119901 ). Therefore, in view of the estimates pre- 
sented in Fig. [12] (see the next Section for details), heating of H2 by X-rays from supernova 
remnants is completely negligible. 



7.2. Heating by shocks 

For galaxies such as NGC3627 and NGC4569, with very high H2 to X-ray luminosity 
ratios, an efficient mechanism has to be invoked to account for the H2 brightness. It has 
been shown recently that galactic shocks can convert a very large fraction of the kineti c 
energy into rotational H2 emission without producing a lot of X-rays (lAppleton et al.ll2006l ). 
We thus propose that large-scale shocks play a major role. NGC3627 is an interacting 
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galaxy in the Leo Triplet, characterized by severe morphological and kinematical distortions 
and, in the center, a massive molecu lar gas concentration in the form of an asymmetric 
bar- like structure (jZhang et al.l Il993l ). with a peak in the stellar velocity dispersion that 
i s shifted by about 3^^ from t he nucleus, and probably strong g as inflow into the nucleus 
(lAfanasiev fc Sirchenkol boosh . NGC4569 is a starburst LINER jMaoz et al.lll998h having 
recently produced a large nu mber of supernova explosions, tri ggering the expari sion of X- 
ray. Ha ( iTschoke et al.ll200ll ) and synchrotron emission lobes (jChyzy et al. 120061). and also 



has a circumnuclear ring of molecular gas with strong non-circular motions (INakanishi et al. 



20051 ). The situation could be similar to NGC4945, where H2 emiss ion is seen to follow th e 
innermost part of a starburst outflow, with an extent of ~ 200 pc (IMoorwood et al.lll996l ). 
These two examples thus support the idea of a dominant excitation by shocks, triggered by 
dynamical perturbations and by a starburst wind, respectively. 

In order to test the possibility that the excess H2 emission be caused by supernova 
remnant shocks, we estimated the required heating efficiencies in the following way. We first 
computed the star formation rate required to account for the excess H2 emission, assum- 
ing constant star formation, if all the mechanical power is converted to the power emitted 
by the sum of the S (0) to S(2) lines (using the Starburst99 population synthesis model of 
Leitherer et al.lll999l ). We then compared this with the star form ation rate estimated from 
the 24 /im luminosity, using the calibration of IWu et al.l (120051 ) for the same initial mass 
function. The ratio of the two rates gives an order-of-magnitude estimate of the H2 heating 
efficiency rj that is needed if supernova shocks are invoked as the dominant heating mecha- 
nism, and is shown in Fig. [T2] as a function of the excess H2 luminosity. In this simplified 
computation, the H2 excess in the star-forming regions B and C of NGC 6822 can easily be 
attributed to supernova remnant shocks, as they would imply efficiencies of at most 3%. 
Besides these regions, supernova remnants are also a sufficient heating source of the excess 
H2 for at least six of the LINER and Sy nuclei. In NGC 5195, in particular, the necessary 
heating efficiency (including the fraction of the total mechanical power injected and ab- 
sorbed in molecular clouds) is very small. Galaxies such as NGC 4450 and NGC 4579, on the 
other hand, unambiguously require much more power than available in supernova remnants. 
The excitation source could be shocks triggered by cloud collisions induced by gravitational 
perturbations, maybe in combination with X-rays. 

NGC 1316 (Fornax A) stands out in Fig. [7] as having the highest Ti temperature (Ti > 
280 K) and the lowest mass fraction in warm phase (< 1%), and the lowest levels up to J=4 
are characterized by a single temperature. In this particular case, H2 may be heated by 
fast shocks caused by the powerful jet that has produced large-scale radio lobes, and that is 
observed in the form of knots in the inner kiloparsec, which may be a sig nature of interaction 
with the interstellar medium, i.e. shocks (iGeldzahler fc Fomalontlll984j ). Excitation by far- 
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ultraviolet radiation and X-rays or by s low shocks would 
at temperatures between 100 and 300 K (IBurton et al.lll99 



i 



roduce a large quari tity of H2 



Maloney et al.lll996l ). which is 



ruled out by the data. NGC 1316 is the only galaxy of the sample for which the excess H2 
may be heated by an AGN jet. 

Of the 16 H2-excess galaxies, all have excitation diagrams consistent with OPi? high t = 3 
(at least in the free-Ti fits), except NGC 1266, NGC 5866 and NGC 4125. NGC 1266 may be 
unusual simply by virtue of its significant optical depth, as discussed in Section 13.31 This 
predominance of OPi? high t = 3, combined with more elevated temperatures than in HII 
nuclei on average (Fig. [7]) would be consistent with the excess H2 emission originating in 
shocks where ortho-para equilibration is fast. It remains however difficult to test this idea, 
and to identify the source of these shocks. 



8. Summary and conclusions 

We present the measurements and results of a survey of the four lowest-energy rotational 
transitions of H2, S(0) to S(3), in a local sample of 57 galaxies, from the SINGS program. 
For three galaxies in this sample, higher-energy transitions, up to S(7), could be measured. 
Characterizing the amount and physical conditions of the warm molecular hydrogen phase 
traced by these lines is of prime interest, because molecular hydrogen represents a major 
mass fraction of the interstellar medium of normal galaxies, and the warm phase itself (gas 
heated to temperatures of ~ 100 to 1000 K) can constitute a substantial fraction of the total 
H2. 

The emission is measured over areas of median size 0.9 kpc, thus including a large 
number of distinct star formation sites and molecular clouds. The sample comprises mostly 
nuclear regions (in 47 massive galaxies and 9 dwarf galaxies), of which 45% are optically 
classified as LINER or Seyfert, as well as 10 extranuclear star-formation complexes within a 
dwarf galaxy and a spiral. With respect to earlier studies of molecular hydrogen emission in 
galaxies, and particularly rotational lines, which had focussed on very bright systems (nearby 
starbursts and AGNs, as well as ultraluminous galaxies), this paper provides results on the 
average properties of warm H2 of relatively faint systems, more representative of the general 
population of galaxies. 

Perhaps the most significant observational results (detailed below) are: (1) the tight 
correlation of the powers emitted by the sum of the S(0) to S(2) lines and by aromatic bands, 
and the fact that the P(S0-S2)/PAH ratio is insensitive to the marked variations in average 
radiation field intensities existing in our sample; (2) the existence of non-equilibrium ortho 
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to para ratios in the rotational levels, that are weakly correlated with the surface brightness 
of the H2 lines. These results call for further modelling in order to be better understood. 

Masses and column densities: The total masses of warm H2 within our apertures 
range between 10^ Mq and close to 3 x 10^ Mq in galaxy nuclear regions. In star formation 
complexes of nearby dwarf galaxies (NGC2915 and NGC6822), we probe warm H2 masses 
down to a few 10^ Mq within equivalent diameters of 60 to 250 pc. The mass densities range 
between 0.2 and 30 Mq pc~^. The column densities that we derive are on average of the same 
order of magnitude as the column densities observed in individual Galactic photodissociation 
regions. For systems in which H2 is predominantly excited by PDRs, assuming they have 
similar characteristics to the Orion Bar, this implies that they fill most of the observing 
beam. 

H2 mass fraction in the warm phase (T > 100 K): Under a conservative assumption 
about the distribution of temperatures, we find that the warm H2 gas makes up between 1% 
and more than 30% of the total H2. For star- forming galaxies, the median mass fraction in 
the warm phase is 10%. The column density nevertheless has a steep inverse dependence 
on the temperature, and we cannot rule out that the unconstrained cool H2 component 
(70 K < T < 100 K) might in some cases dominate the S(0) and S(l) emission and account 
for most of the H2 mass. 

PDR excitation in star-forming regions: In HII nuclei, we observe a remarkably 
narrow range of H2 to aromatic band flux ratio. This result argues for photodissociation re- 
gions providing most of the power used for II2 excitati on, since aromatic bands are known to 



arise predominantly from these regions, as defined by iHoUenbach fc Tielend (119971 ). that in- 
clude all the neutral interstellar medium illuminated by FUV photons. Two main excitation 
mechanisms can be at work simultaneously in PDRs: pumping by FUV photons, followed 
by fluorescent decay to the ground electronic state, and collisions with hydrogen atoms and 
molecules heated by photoelectrons or FUV pumping. Comparison with the predictions of 
PDR models for the ratio of H2 to FUV luminosity indicates a narrow range of average 
physical conditions, Go/n between 0.1 and 1 with the radiation field intensity Gq in units of 
1.6 X 10~'^ergs~^ cm~^ and the hydrogen nucleus density n in units of cm~^. The sum of the 
S(0) to S(2) transitions represents between 2.5 x 10~^ and 7.5 x 10"^ of the total infrared 
power, and on average 30% of the [Sill] line power. The observed temperatures suggest 
that the average Go varies between about 100 and 5000, and the average n between about 
500 and lO^cm"^. This seems to imply that H2 rotational line emission comes mostly from 
molecular clouds illuminated by OB associations. We have seen however that the estimator 
of these parameters is biased (see Section \Q7L\i . and a non- negligible contribution from less 
dense regions with less intense radiation fields is not excluded. 
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Evidence for fluorescence in steir-forniing regions: Previous surveys of rotational 
lines in galaxies had not revealed any departure from thermalization of H2 , which is a conse- 
quence of the relatively low critical densities of the lower rotational transitions. By contrast, 
we find that nearly half the targets in our sample deviate significantly from local ther- 
modynamic equilibrium in having apparent ortho to para ratios OP-RhighT lower than the 
equilibrium value of three. We have seen that this result cannot be an artifact caused by 
extinction effects. Low values of OPi^highT may thus be interpreted as evidence of fluo- 
rescent excitation, which naturally leads to low ortho to para ratios in the excited states, 
occurring in regions of sufficiently low density that ortho-para equilibration by collisions is 
incomplete. The fraction of relatively diffuse molecular gas in normal galaxies could thus be 
far from negligible. In order to test this idea, it would be desirable to obtain independent 
estimates of the average gas density in PDRs. In the absence of any robust estimate of 
this quantity, we used the surface brightness of the sum of the S(0) to S(2) transitions as a 
tracer of the coupled variations of n and Go to show that the data arc compatible with an 
interpretation of low apparent OP-RhighT values in terms of low density PDRs. In particu- 
lar, the lowest values occur preferentially in very late type galaxies, that may have a more 
diffuse interstellar medium than earlier-type galaxies. In the present sample, close to half 
the star- forming targets have low OP-RhighT values. We infer that ffuorescence can be the 
predominant excitation mechanism of rotational H2 lines in normal star-forming galaxies. In 
more active galaxies, however, coUisional excitation in photodissociation regions is likely to 
overtake fluorescence, even though the latter still contributes to gas heating. Alternatively, 
low OPR highT values may be caused by non-equilibrium PDRs, in which initially cold gas, 
recently reached by the photodissociation front and heated, has not had enough time to 
adjust its ortho to para ratio. A disadvantage of this scenario is that it does not explain why 
PDRs would have systematically different properties in targets with OPP highT < 3 and in 
those with OPR highT — 3. We thus emphasize that the cause of non-equilibrium ortho to 
para ratios in the rotational levels is not well understood at present. 

Differences between HII and LINER/Sy nuclei: Despite our large observing beam 
implying that the emission from the immediate vicinity of the nucleus is diluted in the 
emission from extended areas decoupled from nuclear activity, a large fraction of nuclei 
classifled as LINER or Seyfert distinguish themselves from purely star-forming nuclei in 
several ways. In a statistical sense, the temperatures of the warm H2 phase are slightly 
higher, and as a corollary the mass fractions of warm to total H2 are lower (with a median of 
4% instead of 10% for HII nuclei). The correlation between II2 and aromatic band emission 
observed in HII nuclei also breaks down in LINER and Sy nuclei. A large number of them 
have excess emission in the H2 lines, with respect to aromatic bands and to [Sill] , which in 
general is the brightest cooling line in the mid-infrared range (with possible contributions 
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from both PDRs and HII regions as well as X-ray irradiated gas), and, to a lesser degree, 
with respect to the total infrared emission. The fact that less contrast between the different 
nuclear categories is seen in the H2/TIR ratio may partly be due to the fact that estimating 
the 70 fim and 160 /im fluxes within our small apertures requires a large extrapolation. 

Threshold for nonstellar excitation: We propose that nuclear ratios F{SO — S2)/F7.9dust > 
20-1.94^ with observables defined as in our study, are indicative of the LINER and Sy cat- 
egories. It should however be kept in mind that some sources belonging to these classes 
are indistinguishable from HII nuclei. This could be thought of as mainly a distance and 
beam dilution effect; but the H2/PAH ratio behaves contrary to this expectation, since in the 
present sample it shows a correlation with distance, rather than an anti-correlation. None of 
the quantities derived in this paper shows any dependence on distance, with this exception of 
the H2 excess in LINER and Sy nuclei. The fact that it tends to increase with distance could 
be a selection bias, but it also suggests that the H2 excess is in general spatially extended. 
We may also remark that for galaxies like NGC 1377, who se dust emission is interpreted 



as dominated by an extremely young and opaque starburst (IRoussel et al.ll2006l ). the above 
criterion to select LINER and Sy nuclei will not be applicable. 

Shock excitation in LINER/Sy nuclei: We interpret the differences in H2/PAH 
ratios as a genuine excess of H2 emission (and not a deficit of the other tracers while pre- 
serving H2), i.e. as requiring at least one additional mechanism to excite H2 molecules with 
respect to FDR heating. Excitation by nuclear X-rays seems implausible, as models predict 
much lower heating efficiencies than what would be necessary to account for the estimated 
H2 excess, compared with X-ray luminosities derived from Chandra observations. We thus 
favor excess heating by large-scale shocks, caused either by the collective effect of supernovse 
in an aging starburst, or by dynamical perturbations. An order-of-magnitude estimate sug- 
gests that supernova remnant shocks can easily account for the H2 excess of a fraction of 
the LINER and Sy nuclei (for example NGC 5195), but do not provide enough mechanical 
power for the galaxies with the highest H2/PAH ratios. For the latter, shocks triggered by 
dynamical perturbations are the best candidate to supply the excess H2 heating. In one 
case, it is conceivable that the excess H2 emission may be produced by the interaction of a 
nuclear jet with the interstellar medium, namely in NGC 1316 (Fornax A). This target has 
the warmest H2 of the whole sample, with no evidence for a cool (T < 300 K) component, 
the lowest mass fraction in the warm phase, and is remarkable for its large-scale radio lobes. 

Ortho-para thermalization in LINER/Sy nuclei: Consistent with the hypothesis 
that the additional II2 may be caused by shocks in which ortho-para equilibration is fast, 
the excitation diagrams of most galaxies with excess H2 emission are consistent with ortho- 
para thermalization (OP-RhighT = 3), with three exceptions: NGC 1266, where we have 
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seen that the apparent OPi? high t < 3 may be a resuh of high optical depth toward H2, 
NGC5866 and NGC4125. Note that the latter two galaxies are observed edge-on, so that 
OPi? high T < 3 could also be an extinction artefact; this would however imply that the optical 
depth derived from the silicate bands be underestimated by large factors (see Sect. I3.3p . If 
excess H2 emission originated from shocks with OPP high t < 3 in the latter two galaxies, it 
would imply that the gas was initially cold and had not had time to reach equilibrium. 

The rotational H2 lines are most often fainter than the forbidden lines in the same 
wavelength range, in particular [Sill], [Nell], [Nelll] and [SlII] (with no example of an H2- 
dominated line spectrum in the present sample), but are among the dominant coolants of 
molecular gas and provide important constraints on the excitation of the warm molecular 
interstellar medium, where a large mass fraction of the gas resides in normal galaxies. The 
results presented here are assumed to be representative of moderate-luminosity galaxies of 
all types, and can serve as a comparison point for future studies of distant galaxies. 

Support for this work, part of the Spitzer Space Telescope Legacy Science Program, was 
provided by NASA through an award issued by the Jet Propulsion Laboratory, California 
Institute of Technology under NASA contract 1407. We thank Adam Leroy for a useful 
discussion. 
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Table 1. Observed H2 lines. 



transition 


short 


rest A 


spectral 




A 


v=0 


notation 


(/xm) 


order 


(K) 


(10-" s-i) 


J=2-0 


S(0) 


28.219 


LH 14 


510 


2.95 


J=3-l 


S(l) 


17.035 


SH12 


1015 


47.6 


J=4-2 


S(2) 


12.279 


SH17 


1681 


275. 


J=5-3 


S(3) 


9.665 


SLl 


2503 


980. 


J=6-4 


S(4) 


8.025 


SLl 


3473 


2640. 


J=7-5 


S(5) 


6.910 


SL2 


4585 


5880. 


J=8-6 


S(6) 


6.109 


SL2 


5828 


11400. 


J=9-7 


S(7) 


5.511 


SL2 


7196 


20000. 



'^ The rotational upp e r leve l energies were computed from the molecul ar constants given 
bv iHuber fc Herzberg (119791 ) and the transition probabilities are from iBlack fc Dalgarno 
fll976h . 
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Table 2. Targets.^ 



galaxy D class center (J2000) solid angle 

(Mpc) RA DEC (arcsec^) 



N24 


8, 


,2 


dwRrf 


00 


09 


56, 


,31 


-24 


57 


52, 


,0 


294. 


N337 


24, 


.7 


nuc 


00 


59 


49, 


.99 


-07 


34 


42, 


,5 


289. 


N628 


11, 


.4 


nuc 


01 


36 


41, 


.72 


+15 


46 


59, 


,4 


300. 


N855 


9, 


.6 


dwarf 

KX VV CCJ- J- 


02 


14 


03, 


.64 


+27 


52 


40, 


.6 


288. 


N925 


10, 


,1 


IlllC 


02 


27 


17, 


,06 


+33 


34 


43, 


,8 


287. 


N1097 


16, 


.9 


nuc 


02 


46 


18, 


.77 


—30 


16 


30, 


,0 


809. 


N1266 


31, 


.3 


liner 


03 


16 


00, 


.68 


-02 


25 


38, 


,8 


287. 


N1291 


9, 


,7 


liner 


03 


17 


18, 


,50 


-41 


06 


27, 


,7 


287. 


N1316 


26, 


,3 


liner 


03 


22 


41, 


,61 


-37 


12 


28, 


,7 


287. 


N1482 


22, 


.0 


nuc 


03 


54 


38, 


,68 


-20 


30 


08, 


,5 


798. 


N1512 


10, 


.4 


nuc 


04 


03 


53, 


.96 


—43 


20 


55, 


.6 


831. 


N1566 


18, 


,0 




04 


20 


00, 


,53 


-54 


56 


17, 


,4 


306. 


N1705 


5, 


,8 


dwarf 


04 


54 


13, 


.26 


-53 


21 


39, 


,2 


802. 


N2403 


3. 


.5 


nuc 


07 


36 


50, 


.25 


+65 


36 


04, 


.6 


287. 


N2798 


24, 


.7 


nuc 


09 


17 


22, 


.99 


+42 


00 


00, 


.8 


293. 


N2841 


9, 


,8 


liner 


09 


22 


02, 


,75 


+50 


58 


35, 


,5 


291. 


N2915 


2, 


,7 


dwarf 


09 


26 


11, 


.59 


-76 


37 


34, 


,2 


291. 


N2976 


3, 


.5 


dwarf 

KX VV uJ- J- 


09 


47 


15, 


.63 


+67 


55 


00, 


.3 


289. 




Q 
o , 


K 
, o 


sy 




oo 


oo , 






uo 


oo. 


a 
, o 




N3049 


19, 


,6 


nuc 


09 


54 


49, 


,66 


+09 


16 


19, 


,2 


320. 


N3184 


8, 


,6 


nuc 


10 


18 


17, 


.03 


+41 


25 


28, 


.3 


287. 


N3190 


17, 


.4 


liner 


10 


18 


05, 


.71 


+21 


49 


56, 


.8 


304. 


N3198 


9, 


,8 


nuc 


10 


19 


55, 


,10 


+45 


32 


59, 


,6 


287. 


N3265 


20, 


.0 


nuc 


10 


31 


06, 


.85 


+28 


47 


49, 


.1 


301. 


Mrk33 


21, 


.7 


dwarf 


10 


32 


32, 


.02 


+54 


24 


05, 


.4 


315. 


N3351 


9. 


.3 


nuc 


10 


43 


57, 


.83 


+11 


42 


10, 


.9 


761. 


N3521 


9, 


,0 


liner 


11 


05 


48, 


,51 


-00 


02 


10, 


,1 


285. 


N3621 


6, 


,2 


liner 


11 


18 


16, 


.62 


-32 


48 


49, 


,9 


291. 


N3627 


8, 


.9 


sy 


11 


20 


15, 


.04 


+12 


59 


31, 


.0 


301. 


N3773 


12, 


.9 


dwarf 


11 


38 


12, 


,91 


+12 


06 


43, 


,0 


285. 


N3938 


12, 


,2 


nuc 


11 


52 


49, 


,45 


+44 


07 


15, 


,0 


292. 


N4125 


21, 


.4 


liner 


12 


08 


05, 


.73 


+65 


10 


27, 


.9 


294. 


N4254 


20, 


.0 


nuc 


12 


18 


49, 


.71 


+14 


25 


01, 


.4 


299. 


N4321 


20, 


,0 


nuc 


12 


22 


55, 


.08 


+15 


49 


18, 


,2 


795. 


N4450 


20, 


.0 


liner 


12 


28 


29, 


.51 


+17 


05 


04, 


.9 


285. 
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Table 2 — Continued 



galaxy D class center (J2000) solid angle 

(Mpc) RA DEC (arcsec^) 



N4536 


25.0 


nuc 


12 


34 


27 


22 


+02 


11 


14 


3 


777. 


N4552 


20.0 


liner 


12 


35 


39 


71 


+ 12 


33 


21 


9 


285. 


N4559 


11.6 


nuc 


12 


35 


57 


79 


+27 


57 


36 


7 


295. 


N4569 


20.0 


liner 


12 


36 


49 


91 


+13 


09 


46 


8 


287. 


N4579 


20.0 


sy 


12 


37 


43 


70 


+ 11 


49 


07 


1 


295. 


N4594 


13.7 


liner 


12 


39 


59 


58 


-11 


37 


22 


3 


289. 


N4625 


9.5 


dwarf 


12 


41 


52 


65 


+41 


16 


26 


5 


294. 


N4631 


9.0 


nuc 


12 


42 


07 


84 


+32 


32 


34 


8 


282. 


N4725 


17.1 


sy 


12 


50 


26 


71 


+25 


30 


03 


2 


294. 


N4736 


5.3 


liner 


12 


50 


53 


23 


+41 


07 


13 


6 


298. 


N4826 


5.6 


liner 


12 


56 


43 


79 


+21 


41 


00 


7 


292. 


N5033 


13.3 


sy 


13 


13 


27 


67 


+36 


35 


38 


2 


311. 


N5055 


8.2 


liner 


13 


15 


49 


50 


+42 


01 


46 


6 


313. 


N5194 


8.2 


sy 


13 


29 


52 


94 


+47 


11 


44 


2 


330. 


N5194jiL'= 


8.2 


hii 


13 


29 


49 


54 


+47 


13 


28 


6 


256. 


"IN. T c 1 n /I 13 b 


Q O 

o.z 


nil 




oU 




oy 


+4/ 




OU 


n 

y 


071 


N5194.C^ 


8.2 


hii 


13 


30 


00 


02 


+47 


11 


11 


9 


352. 


N5194J3'^ 


8.2 


hii 


13 


30 


02 


76 


+47 


09 


53 





338. 


N5194J]'^ 


8.2 


hii 


13 


29 


56 


84 


+47 


10 


45 


9 


234. 


N5194J'^ 


8.2 


hii 


13 


29 


52 


93 


+47 


12 


38 


7 


297. 


N5194.G'^ 


8.2 


hii 


13 


29 


43 


99 


+47 


10 


20 


7 


285. 


N5195 


8.2 


liner 


13 


29 


59 


80 


+47 


16 


00 


1 


323. 


N5713 


26.6 


nuc 


14 


40 


11 


57 


-00 


17 


19 


1 


292. 


N5866 


12.5 


liner 


15 


06 


29 


58 


+55 


45 


46 


7 


286. 


N6822^'= 


0.6 


dwarf 


19 


44 


52 


84 


-14 


43 


09 


8 


308. 


N6822.B'= 


0.6 


dwarf 


19 


44 


50 


57 


-14 


52 


49 





288. 


N6822.C'^ 


0.6 


dwarf 


19 


44 


48 


67 


-14 


52 


26 


5 


304. 


N6946 


5.5 


nuc 


20 


34 


51 


79 


+60 


09 


12 


6 


803. 


N7331 


15.7 


liner 


22 


37 


04 


11 


+34 


24 


57 


5 


263. 


N7552 


22.3 


nuc 


23 


16 


10 


56 


-42 


35 


04 





722. 


N7793 


3.2 


dwarf 


23 


57 


49 


79 


-32 


35 


28 


4 


296. 



^Distances are from Kennicutt at al.l (j2003l ). Nuclear classifications derived from optical spectroscopy, 
published by ISmith et al.l |2007bl ). were modified for the following targets. N1097, N1512, N4321 and 
N7552: Since our aperture includes a bright star-forming ring, which dominates the dust and line emission, 
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we adopt the HII class instead o f the LINER cl ass. N2841, N4552, N4569, N4594, N4826 and N5195: The 
Sy class w as changed to LINER |Ho et al.lll997l) . N3198 and N3938: The LINER class was changed to HII 
( Ho et al.lll997.) . The dwarf galaxy class is here arbitrarily defined by a total stellar mass, estimated as by 



Lee et al. 



( 20061) . below IO^ '^Mg 



'^A total of 11 extranuclear regions were observed in high-resolution spectroscopy in NGC 5194, from which 
we included 7 in our sample. 

'^Several locations were observed in high-resolution spectroscopy in this galaxy. We present here results 
for three amo ng the brightest regions at 7.9 and 24 /xm. N6822_A corre sponds to the HII region Hubble V 
|Hubblelll925[) . and N6822_C to the HII region K7 jxinman et al.lll979[) . 



Table 3. Fluxes.^ 



galaxy 


S(0) 


S(l) 


S(2) 


S(3) 


[Sill] 


Fy.sdust 


F24 


TIR 








(10-18 Wm-2 


) 




(Jy) 


(Jy) 


(10-15 Wm-2) 


N24 


6.6 ±2.6 


2.2 ±2.4 


3.2 ±2.0 


10.3 ±6.4 


36.7 ±3.6 


0.010 


0.010 


18. ±3. 


N337 


11.6 ±3.8 


21.0 ±4.9 


15.4 ±9.9 


9.7 ±5.2 


186.1 ±7.0 


0.048 


0.087 


89. ± 9. 


N628 


7.9 ±2.5 


13.7 ±5.3 


0.0 ±3.0 


5.7 ±8.6 


47.8 ± 2.8 


0.026 


0.026 


39. ± 3. 


N855 


2.9 ±0.9 


11.2±6.5 


0.0±5.1 


22.4 ±7.2 


70.2 ±3.8 


0.017 


0.039 


53. ± 7. 


N925 


3.8 ± 1.7 


7.8 ±2.9 


3.7 ±2.1 


0.0 ±8.8 


81.0 ±4.5 


0.021 


0.022 


37. ±4. 


N1097 


213.1 ±42.6 


726.1 ±43.3 


293.6 ± 29.4 


423.0 ±23.4 


4921.1 ± 121.9 


1.139 


3.763 


2214. ± 227. 


N1266'' 


16.4 ±8.5 


148.5 ±6.6 


121.8 ±7.1 


189.8 ± 11.6 


153.5 ±41.5 


0.048 


0.579 


438. ± 30. 


N1291 


3.7 ± 1.8 


29.8 ±4.8 


13.9 ±8.2 


31.6 ± 10.0 


75.1 ±5.8 


0.013 


0.049 


60. ± 7. 


N1316 


1.5 ±0.8 


36.0 ±6.1 


18.4 ±5.9 


84.0 ±9.9 


102.4 ±6.4 


0.025 


0.094 


105. ±11. 


N1482 


106.8 ±35.1 


424.0 ± 53.4 


184.0 ± 22.8 


207.5 ± 26.9 


3954.0 ± 75.1 


1.163 


2.974 


1752. ± 154. 


N1512 


25.9 ±2.9 


83.6 ±8.9 


27.2 ±10.4 


48.3 ± 13.6 


390.7 ± 13.9 


0.091 


0.150 


203. ± 24. 


N1566 


24.0 ±2.2 


129.5 ±7.1 


55.3 ±4.8 


91.6 ± 17.6 


101.9 ±7.8 


0.098 


0.237 


235. ±19. 


N1705 


5.6 ±3.9 


2.2 ± 1.3 


0.0 ± 10.3 


11.6 ±4.0 


86.3 ± 12.0 


0.007 


0.028 


40. ±5. 


N2403 


5.5 ± 1.5 


6.6 ±3.1 


2.0 ±2.0 


4.6 ±8.2 


82.0 ±5.9 


0.029 


0.022 


35. ± 2. 


N2798 


43.2 ±25.3 


207.8 ±22.1 


90.7 ±9.2 


105.2 ± 12.1 


840.4 ±42.2 


0.354 


1.723 


936. ±67. 


N2841 


3.1 ±1.4 


9.6 ±2.5 


4.1 ±4.5 


28.0 ± 10.8 


81.1 ±5.5 


0.010 


0.027 


48. ±4. 


N2915 


0.8 ±0.6 


1.8 ±2.1 


3.8 ± 1.8 


0.0 ±8.0 


40.1 ±3.7 


0.005 


0.019 


25. ± 3. 


N2976 


8.9 ±2.4 


15.7±3.0 


4.9 ±3.5 


5.2 ±4.2 


66.1 ±5.4 


0.026 


0.055 


77. ± 6. 


N3031 


6.3 ±3.6 


46.1 ±6.1 


26.1 ±8.1 


65.3 ±8.0 


217.9 ±9.7 


0.081 


0.337 


234. ± 17. 


N3049 


6.4 ±2.3 


23.7 ±3.8 


11.7±6.5 


11.1 ±5.2 


197.1 ± 11.9 


0.051 


0.258 


141. ±13. 


N3184 


9.8 ±2.3 


22.6 ±2.9 


6.3 ±3.5 


0.0 ± 15.9 


105.0 ±4.2 


0.031 


0.082 


103. ± 13. 


N3190 


18.1 ±2.9 


75.3 ±8.0 


20.9 ±6.4 


71.6 ± 12.6 


67.1 ±7.3 


0.051 


0.068 


137. ±23. 


N3198 


13.2 ±4.9 


32.1 ±5.4 


10.2 ±3.5 


18.7 ±8.5 


96.5 ±6.6 


0.059 


0.260 


200. ± 20. 


N3265 


9.3 ±3.2 


30.6 ±3.5 


11.4 ± 7.9 


15.3 ±7.4 


159.1 ±7.7 


0.062 


0.193 


119. ±9. 


Mrk33 


13.9 ±5.6 


29.3 ±3.7 


11.5 ±5.3 


45.7 ±8.6 


245.4 ± 10.9 


0.076 


0.583 


239. ± 13. 


N3351 


62.8 ± 10.5 


218.4 ±20.6 


88.3 ± 10.1 


169.8 ± 19.9 


1858.9 ±35.7 


0.366 


1.385 


917. ± 93. 



Table 3 — Continued 



galaxy 


S(0) 


S(l) S(2) 


S(3) 


[Sill] Fr.Qdust F24 










(Jy) (Jy) (10- 



TIR 

i5Wm-2) 



N3521 


1 7 C _|_ Q C 
-L 1 .0 —1— O.'J 




12.3 ± 4.0 


ZjU . 'J —1— 0.0 


1 c;7 c. -(- 4 « 


101 


1 06 


129. ± 10, 


Ly OU^ J. 




00.0 ZC O.Zi 


171 -1-4 7 




±<JVJ.O ZIZ 1 .(J 


07^ 

U.W 1 


\j.\j\jO 


00. zn u. 


N3627 


31 2 + 3 Q 


tJ-LO<W —1— J. ^ • ^ 


141.7 ± 5.2 


onq 1-1-12 6 


162 1 ±93 


0.190 


0.406 


518. ± 31, 


N3773 


4.8 ± 3.4 


7.3 ± 3.2 


3.7 ± 2.5 


8.3 ± 8.6 




0.024 


0.087 


67. ± 7, 


N3q38 


8.0 ± 1.2 


13.4 ± 3.3 


4.2 ± 3.2 


0.0 ± 8.0 


57.7 ± 4.2 


0.027 


0.030 


51. ± 5. 


N4125 


2.6 ± 1.6 


17.6 ± 4.7 


12.6 ± 6.4 


18 ± 8 8 


46 8 ± 5 9 


0.005 


0.015 


21. ± 3. 


N4254 


17.6 ± 1.3 


87 3 ± 9 5 


46 6 ± 9 2 


32 5 ± 8 8 


409 ± 10 6 


0.160 


0.203 


231 ± 23 


N4321 


78.2 ± 11.7 




1 Qo 1 -1- 24 6 


1 C4 5 -1- 1 c 


1 507 Q -(- 07 2 


0.392 


0.730 


778. ± 90, 


N4450 


10.3 ± 1.6 


91.4 ± 5.2 




oq n -(- 1 1 


00 fi -1- q 1 


0.011 


0.028 


55. ± 9. 


N4536 


1 no 7 _|_ 00 q 


41 1 c _|_ 4n c 


1 7c c -1 o-? Q 

-L f t_J.\J _1_ .'.J 


216 8 ± 20 7 


3207 4 ± 70 


0.813 


2.401 


1 4f)f) -1- 1 1 s 

t _l_ J. J. (Jf 1 


N4552 


1.8 lb 1.4 


1.1 ± 1.5 


0.0 =b 3.9 


9.1 ± 4.6 


13.4 ± 3.9 


0.004 


0.017 


6. ± 1. 


N4559 


13.6 ± 2.1 


17.9 ± 3.5 


4.0 ± 1.5 


OQ 4 -L 1 Q 7 


81.1 ± 3.5 


0.034 


0.030 


47. ± 5, 




41.0 =t D.y 


OiO.i ± 0. ( 


1 c;o ./] 1 c: 
lOZ.4 zt O.Z 


ouy.4 zt io. ( 


0«a r; 4- 1 7 

ZOO.O It i ( .0 


1 oa 


n A 7P^ 

U.4 ( 


"XPrK 4- /in 

oDO. ± 4U, 


N4579'* 


8.9 ±2.3 


164.1 ±5.2 


105.8 ± 12.4 


256.3 ± 10.0 


175.3 ± 12.2 


0.032 


0.135 


160. ± 18, 


N4594 


1.9 ±2.3 


6.4 ±3.8 


4.2 ±2.7 


4.6 ±4.0 


126.1 ±7.3 


0.017 


0.068 


70. ± 8, 


N4625 


8.0 ± 1.6 


9.6 ± 1.9 


5.4 ±4.3 


5.0 ± 11.4 


67.5 ±4.8 


0.024 


0.024 


33. ± 5. 


N4631 


58.2 ±6.4 


123.6 ± 16.1 


60.4 ±5.6 


46.8 ±9.1 


1350.4 ±22.7 


0.338 


0.661 


636. ±61, 


N4725 


11.7± 1.7 


37.1 ±5.1 


19.9 ±8.1 


38.0 ±9.0 


24.9 ±3.2 


0.011 


0.035 


56. ±9. 


N4736 


33.2 ±6.8 


250.8 ± 16.4 


100.1 ±10.4 


221. 7± 13.6 


219.5 ±21.4 


0.293 


0.636 


678. ± 50, 


N4826 


71.3 ±6.5 


345.3 ± 15.8 


150.9 ± 10.9 


210.6 ± 10.3 


927.6 ± 15.0 


0.343 


0.612 


756. ± 92, 


N5033 


36.6 ±3.5 


182.0 ± 10.4 


63.5 ±3.1 


126.9 ± 19.1 


530.9 ± 10.0 


0.177 


0.239 


311. ±45, 


N5055 


44.0 ±2.4 


158.0 ±9.1 


52.0 ±6.7 


80.2 ± 10.6 


295.2 ±9.0 


0.148 


0.166 


300. ±37, 


N5194 


17.7 ±2.8 


134.4 ±9.1 


75.7 ±4.7 


187.3 ± 17.7 


634.9 ± 12.2 


0.168 


0.291 


351. ±39, 


N5194^ 


5.8 ± 1.3 


16.4 ±2.7 


4.3 ±2.6 


5.4 ± 1.8 


27.8 ±3.0 


0.030 


0.036 


53. ±2. 


N5194_B 


18.1 ±2.2 


43.1 ±8.2 


12.0 ±4.3 


18.1 ± 10.2 


339.7 ±8.5 


0.094 


0.193 


201. ±13, 


N5194.C 


23.2 ±2.2 


43.4 ±6.7 


5.7 ±4.5 


14.1 ±4.6 


223.6 ±6.8 


0.079 


0.115 


155. ± 7. 
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galaxy 


S(0) 


S(l) 


S(2) 


S(3) 


[Sill] 


Fy.gdust 


F24 


TIR 








(10-18 Wm-2) 






(Jy) 


(Jy) 


(10-15 Wm-2) 


N5194_D 


14.7 ±2.8 


30.3 ±6.5 


3.9 ±3.1 


8.4 ±3.0 


222.9 ±6.2 


0.058 


0.118 


130. ± 14. 


N5194_E 


22.7± 1.7 


47.9 ±6.0 


19.2 ±7.1 


23.3 ±3.1 


222.2 ±4.8 


0.072 


0.089 


116. ±8. 


N5194_F 


19.4±2.2 


39.2 ±5.2 


15.4 ± 9.1 


21.6 ± 14.0 


271.7 ±5.8 


0.075 


0.110 


132. ±4. 


N5194_G 


22.2 ±3.0 


49.6 ±5.4 


23.6 ±2.7 


27.7± 11.2 


425.6 ±8.5 


0.093 


0.192 


173. ± 15. 


N5195 


54.2 ±20.9 


310.4± 12.3 


127.7±6.5 


275.0 ± 17.2 


127.8 ± 19.6 


0.271 


0.835 


600. ± 55. 


N5713 


27.4 ± 3.6 


151.6 ± 14.3 


50.9 ± 7.9 


85.6 ± 7.0 


798.7 ± 14.8 


0.268 


0.738 


500. ± 54. 


N5866 


15.5± 1.5 


90.0 ±5.5 


39.1 ±7.0 


37.2 ±6.5 


78.8 ±3.9 


0.038 


0.052 


169. ±26. 


N6822_A 


8.1 ±4.3 


21.1 ±4.3 


16.0 ±4.1 


25.9 ±3.8 


187.1 ± 11.2 


0.036 


0.517 


283. ± 24. 


N6822_B 


7.8 ± 1.3 


11.2±3.6 


8.8 ±2.8 


10.9 ±5.1 


19.0 ±2.8 


0.010 


0.037 


44. ± 5. 


N6822_C 


4.4 ± 1.9 


8.0 ±5.2 


6.2 ±5.4 


8.0 ±2.9 


16.0 ±3.6 


0.009 


0.100 


76. ± 7. 


N6946 


227.7±61.9 


636.9 ±33.5 


271.6 ± 17.1 


286.3 ±29.3 


3203.9 ±91.9 


1.314 


4.234 


2354. ± 193. 


N7331 


16.4 ± 1.6 


49.6 ±5.0 


13.8 ± 1.9 


23.4 ±8.3 


213.9 ±4.8 


0.074 


0.091 


125. ± 10. 


N7552 


257.8 ± 106.9 


565.3 ±42.6 


245.9 ±22.5 


315.9 ± 14.8 


3619.9 ± 165.2 


1.495 


7.840 


3237. ±219. 


N7793 


7.0 ± 1.7 


15.0 ±3.2 


11.7 ±4.6 


5.3 ±3.7 


79.7 ±4.3 


0.041 


0.046 


71. ±5. 



''See Sections O and For NGC 1266, NGC 4569 and NGC 4579, it was possible to measure higher-level 
transitions of H2, which are provided in the following notes. 

''The S(4) to S(7) line fluxes are respectively: (102.7 ± 17.4), (241.9 ± 19.9), > 18.5 and (192.2 ± 28.6) x 
10-1* Wm-^. The S(6) transition was not accurately measurable because it was observed at low spectral 
resolution on the blue shoulder of a bright 6.2 /im aromatic band. 

'=TheS(4) toS(7) line fluxes are respectively: (52.6±27.7), < 252.6, < 61.3 and (50.3±24.2) x lO-i* Wm-^. 
The S(5) and S(6) transitions are upper limits because S(5) is contaminated by the [Aril] line at 6.985 /im 
and S(6) is diluted in the blue shoulder of the 6.2 /im aromatic band. 



'^Tho S(4) to S(7) line fluxes are respectively: (48.7 ± 10.7), < 324.9, < 147.2 and (83.3 ± 37.5) x 
10~^® Wm~^. The same remark as for NGC4569 applies. 



I 

CO 

I 
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Table 4. Derived temperatures, column densities and O-P-RhighT-^ 



galaxy 




/i 


T2 


/2 OPRhishT iVtot(T>Ti) Ntot{OPR = 3) 




(K) 




(K) 


(lO^o mol. cm-2) /N^otiOPR < 3) 



N24 


97. 


n 99941 


769. 


00059 


1 40 + n 01 


4.232 


1.75 




78 


n 99980 


675 


00020 


1 Qf) + 57 


15.210 


1.46 




Qfl 


n QQQ71 


888 

000. 


00029 


1 4c _|_ n 07 


7 314 


1 73 

-L . 1 


N337 


Ififl 

xxJxJ. 


9779*^ 


402. 


02207 


1 (12 + 1 5 


1 745 

-L • 1 rkO 


2.13 




128. 


n 97803 


352 


021 97 


1 01 _|_ n OK 


3 218 


97 






9871 1 


464. 


01 289 


n cf; _(_ f) 14 


1.824 


2.19 


N628 


1 1Q 

J. J. C?. 


99fi78 


3Q3 


00322 





2 598 






146. 


n 987^i0 


578 


01 250 





508 




N925 


123. 




494. 


f)0527 


3 


1.050 






99. 


99648 


421. 


00352 


3 


2.780 




N1097 


157. 


98027 


422. 


01973 


no 1 n 

^.\JtJ _1_ \J.(J^J 


11.837 


1.29 






u. y yuoo 




U.UUc/U < 


Q 


fi2 Q2'i 




N1266'' 


247. 


83691 


400. 


15619 

\_/ • -L <J U -L t/ 


1 70 + n 29 


1.231 


1.19 




204. 


82'i88 


400 


16935 


1 an + f) 1 8 


1 594 


84 


N1291 


184 


9'S420 


484 


04580 


3 


41 5 






±4: ± . 


U. iJOyJ i. yJ 


41 4 











N1316 


300. 


96436 


1456. 


03564 


3 


0.114 






> 278. 


> 93750 


> 900. 


< 06250 


3 


0.121 




Ml /too 


1 71 


u.yDy4o 




U.UoUo / 


1.00 ± U.ZO 




i.OO 




137. 


0.96699 


357. 


0.03301 


2.12 ±0.43 


8.848 


1.49 




175. 


0.98252 


469. 


0.01748 


1.58 ±0.25 


5.378 


1.40 


N1512 


138. 


0.98887 


415. 


0.01113 


3 


1.883 






118. 


0.98955 


375. 


0.01045 


3 


3.338 




N1566 


177. 


0.96436 


441. 


0.03564 


2.24 ±0.37 


2.749 


1.19 




115. 


0.97305 


359. 


0.02695 


3 


8.517 




N1705 


127. 


0.99453 


811. 


0.00547 


0.49 ± 0.29 


0.494 


3.78 




92. 


0.99844 


711. 


0.00156 


0.65 ±0.14 


2.050 


0.95 




114. 


0.99727 


936. 


0.00273 


0.49 ±0.11 


0.876 


2.33 


N2403 


109. 


0.99844 


449. 


0.00156 


3 


2.486 






99. 


0.99873 


419. 


0.00127 


3 


4.138 




N2798 


179. 


0.96182 


409. 


0.03818 


1.92 ±0.23 


5.169 


1.32 




144. 


0.95518 


359. 


0.04482 


2.17 ±0.43 


8.460 


1.16 




184. 


0.97812 


472. 


0.02188 


1.67 ±0.27 


5.419 


1.27 


N2841 


137. 


0.99678 


1030. 


0.00322 


3 


0.709 






132. 


0.99707 


970. 


0.00293 


3 


0.844 
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Table 4 — Continued 



galaxy Ti 


/i T2 


/2 OPRi^ighT iVtot(r>Ti) Mot(OPi? = 3) 


(K) 


(K) 


(lO^o mol. cm-2) /Ntot{OPR < 3) 



N2915 


169. 


0, 


.97090 


763. 


0.02910 


1.05 ±0.47 


0, 


,096 


1.74 




116. 


0, 


.98838 


669. 


0.01162 


1.30 ±0.28 


0, 


,305 


0.55 




151. 


0, 


.98350 


880. 


0.01650 


1.00 ±0.20 


0, 


,150 


1.13 


N2976 


134. 





,99199 


376. 


0.00801 


1.84 ±0.30 


2, 


,095 


1.43 




86. 


0, 


.99619 


300. 


0.00381 


3 


12, 


,046 




N3031 


178. 


0, 


.95908 


523. 


0.04092 


3 


0, 


,700 






126. 


0, 


.96387 


433. 


0.03613 


3 


1, 


,699 




N3049 


174. 


0, 


,96611 


403. 


0.03389 


1.60 ±0.20 


0, 


,740 


1.50 




141. 


0, 


.96113 


353. 


0.03887 


1.83 ±0.36 


1, 


,238 


1.17 




179. 


0, 


.98076 


465. 


0.01924 


1.39 ±0.22 


0, 


,773 


1.55 


N3184 


127. 


0, 


.99502 


440. 


0.00498 


3 


2, 


,631 






118. 


0, 


,99648 


449. 


0.00352 


3 


3 


,805 




N3190 


148. 


0, 


,99355 


574. 


0.00645 


3 


3, 


,286 






142. 


0, 


,99385 


554. 


0.00615 


3 


3 


,883 




N3198 


128. 


0, 


,99316 


415. 


0.00684 


3 


3 


,438 






111. 


0, 


,99385 


375. 


0.00615 


3 


6, 


,283 




N3265 


156. 


0, 


,98105 


408. 


0.01895 


2.02 ±0.35 


1, 


,421 


1.30 




95. 


0, 


,98975 


325. 


0.01025 


3 


7, 


,006 




Mrk33 


124. 





,99736 


624. 


00264 

\J . \J\J ZjKJ^ 


3 


3 


,844 






115. 


0, 


,99775 


584. 


0.00225 


3 


5, 


,478 




N3351 


141. 


0, 


,98633 


437. 


0.01367 


3 


4, 


,575 






114. 


0. 


,98867 


387. 


0.01133 


3 


9, 


,893 




N3521 


118. 


0, 


,99531 


414. 


0.00469 


3 


5, 


,715 






91. 


0, 


,99687 


354. 


0.00313 


3 


18, 


,837 




N3621 


126. 


0, 


,99541 


532. 


0.00459 


3 


4, 


,490 






110. 


0, 


,99609 


472. 


0.00391 


3 


8, 


,195 




N3627 


220. 


0, 


,91836 


433. 


0.08164 


2.42 ±0.19 


2, 


,672 


1.12 




124. 


0, 


,90469 


344. 


0.09531 


3 


7, 


,324 




N3773 


115. 


0, 


,99678 


480. 


0.00322 


3 


1, 


,671 






90. 


0, 


,99805 


410. 


0.00195 


3 


5, 


,446 




N3938 


118. 


0, 


,99629 


425. 


0.00371 


3 


2, 


,603 






102. 


0, 


,99697 


377. 


0.00303 


3 


5, 


,116 




N4125 


211. 


0, 


,92803 


481. 


0.07197 


1.75 ±0.24 


0, 


,238 


1.35 




163. 


0, 


,91826 


422. 


0.08174 


1.90 ±0.37 


0, 


,377 


1.15 




212. 


0, 


,95625 


555. 


0.04375 


1.56 ±0.25 


0, 


,262 


1.35 
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Table 4 — Continued 



galaxy Ti 


/i 


T2 


/2 OPRhishT A^tot(T>ri) iVtot(0Pi? = 3) 


(K) 




(K) 


(1020 mol. cm-2) /iVtot(OP-R < 3) 



N4254 


200. 


0.94111 


384. 


0.05889 


1.49±0.11 


1.763 


1.53 




167. 


0.90293 


337. 


0.09707 


1.59 ±0.29 


2.368 


1.14 




206. 


0.96982 


444. 


0.03018 


1.36 ±0.22 


1.812 


1.45 


N4321 


164. 


0.97480 


427. 


0.02520 


1.76 ±0.29 


4.001 


1.41 




131. 


0.97627 


374. 


0.02373 


2.07 ±0.43 


7.511 


1.27 




168. 


0.98477 


493. 


0.01523 


1.49 ±0.25 


4.258 


1.46 


N4450 


190. 


0.96992 


536. 


0.03008 


3 


1.178 






172. 


0.96152 


476. 


0.03848 


3 


1.467 




N4536 


166. 


0.97344 


413. 


0.02656 


1.90 ±0.28 


5.602 


1.35 




133. 


0.97334 


362. 


0.02666 


2.23 ± 0.45 


10.170 


1.35 




170. 


0.98437 


476. 


0.01563 


1.62 ±0.27 


5.908 


1.40 


N4552 


126. 


0.99463 


986. 


0.00537 


0.84 ±0.56 


0.449 


2.61 




89. 


0.99873 


865. 


0.00127 


1.09 ±0.24 


2.156 


0.54 




108. 


0.99766 


1138. 


0.00234 


0.87 ±0.19 


0.959 


1.30 


N4559 


111. 


0.99951 


775. 


0.00049 


3 


6.202 






109. 


0.99961 


825. 


0.00039 


3 


6.967 




N4569'= 


184. 


0.91348 


400. 


0.08564 


3 


4.062 






138. 


0.95226 


400. 


0.04755 


3 


9.055 




N4579d 


286. 


554*^9 


400. 


4sqqq 


on 1 ri 

^,Kj\J _1_ L/.OlJ 


0.534 


1.00 




68. 


0.98551 


400. 


0.01439 


3 


21.592 




N4594 


175. 


0.96621 


441. 


0.03379 


1.49 ±0.23 


0.232 


1.57 




138. 


0.96709 


387. 


0.03291 


1.71 ±0.35 


0.421 


1.06 




178. 


0.97969 


509. 


0.02031 


1.28 ±0.21 


0.251 


1.64 


N4625 


133. 


0.99268 


411. 


0.00732 


1.27 ±0.28 


1.856 


1.85 




107. 


0.99453 


361. 


0.00547 


1.65 ±0.35 


4.144 


1.20 




136. 


0.99551 


475. 


0.00449 


1.05 ±0.20 


1.964 


2.38 


N4631 


154. 


0.98174 


386. 


0.01826 


1.36 ±0.19 


9.852 


1.75 




125. 


0.98174 


338. 


0.01826 


1.64 ±0.32 


18.343 


0.96 




159. 


0.98955 


445. 


0.01045 


1.14 ±0.19 


10.160 


2.01 


N4725 


155. 


0.98223 


481. 


0.01777 


1.97 ±0.50 


1.814 


1.31 




88. 


0.99482 


381. 


0.00518 


3 


13.562 




N4736 


180. 


0.96455 


478. 


0.03545 


3 


3.822 






151. 


0.95732 


418. 


0.04268 


3 


5.791 




N4826 


176. 


0.96523 


426. 


0.03477 


2.05 ±0.30 


8.766 


1.27 




72. 


0.99570 


328. 


0.00430 


3 


223.50 
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Table 4 — Continued 



galaxy 


(K) 




/l 


T2 

(K) 




/2 


OPi? high T 


NtotiT > Ti) 
(10^0 mol. cm-2) 


^tot 

/A^tot 


{OPR = 3) 
[OPR < 3) 


N5033 


157. 


0, 


.97900 


438. 


0, 


.02100 


3 


5.255 








134. 


0, 


.97920 


398. 


0, 


,02080 


3 


8.435 






N5055 


142. 


0, 


.98242 


386. 


0, 


,01758 


3 


7.586 








114. 


0, 


.98506 


346. 


0, 


.01494 


3 


16.336 






N5194 


180. 


0, 


.95732 


521. 


0, 


.04268 


3 


1.753 








127. 


0, 


,96123 


431. 


0, 


,03877 


3 


4.171 






N5194_A 


133. 


0. 


.98672 


349. 


0, 


,01328 


3 


1.469 








111. 


0, 


.98848 


319. 


0, 


.01152 


3 


2.951 






N51943 


128. 


0, 


.99150 


373. 


0, 


,00850 


3 


4.933 








109. 


0, 


.99287 


343. 


0, 


,00713 


3 


9.421 






N5194.C 


120. 


0, 


.99873 


458. 


0, 


,00127 


3 


6.824 








118. 


0, 


.99873 


448. 


0, 


.00127 


3 


7.448 






N5194J) 


123. 


0, 


.99824 


429. 


0, 


,00176 


3 


4.144 








121. 


0, 


.99824 


419. 


0, 


,00176 


3 


4.572 






N5194_E 


142. 


0, 


.98896 


406. 


0, 


.01104 


1.79 ±0.33 


5.653 




1.44 




63. 


0, 


.99922 


313. 


0, 


.00078 


3 


235.5 






N5194J 


137. 


0, 


,99102 


417. 


0, 


,00898 


1.91 ±0.40 


4.127 




1.37 




85. 


0, 


.99648 


331. 


0, 


,00352 


3 


26.208 






N5194.G 


148. 


0, 


,98594 


418. 


0, 


.01406 


1.64 ±0.31 


4.041 




1.53 




118. 


0, 


,98857 


367. 


0, 


,01143 


2.03 ±0.43 


8.457 




1.21 




151. 


0, 


,99160 


483. 


0, 


,00840 


1.37 ±0.24 


4.337 




1.74 


N5195 


164. 


0. 


,97441 


468. 


0, 


,02559 


3 


6.769 








134. 


0, 


,97305 


408. 


0, 


.02695 


3 


11.899 






N5713 


162. 


0, 


,97266 


413. 


0, 


,02734 


3 


3.901 








133. 


0, 


,96875 


363. 


0, 


,03125 


3 


6.684 






N5866 


194. 


0, 


,94678 


393. 


0, 


,05322 


1.88 ±0.15 


1.676 




1.33 




160. 


0, 


,92070 


344. 


0, 


.07930 


2.03 ±0.37 


2.379 




1.10 




200. 


0, 


,97158 


453. 


0, 


,02842 


1.68 ±0.27 


1.731 




1.25 


N6822_A 


158. 


0, 


,98027 


504. 


0, 


,01973 


1.53 ±0.42 


1.094 




1.56 




120. 


0, 


,98721 


442. 


0, 


,01279 


1.89 ±0.41 


2.561 




1.31 




156. 


0, 


,98799 


581. 


0, 


.01201 


1.33 ±0.25 


1.306 




1.76 


N6822_B 


140. 


0, 


,99043 


476. 


0, 


,00957 


1.29 ±0.35 


1.609 




1.82 




109. 


0, 


,99404 


417. 


0, 


,00596 


1.63 ±0.36 


3.946 




1.30 




139. 


0, 


,99414 


549. 


0, 


.00586 


1.09 ±0.21 


1.867 




2.00 


N6822.C 


148. 


0, 


,98604 


478. 


0, 


,01396 


1.33 ±0.35 


0.716 




1.77 
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Table 4 — Continued 





Ti 
(K) 




T2 
(K) 




1 -it high i 


^ * tot 1 ) 

(10^0 mol. cm-2) 


-^'totv^-^ — ^/ 

/Ntot{0PR<3) 




115. 


0.99092 


419. 


0.00908 


1.67 ±0.34 


1.690 


1.04 




148. 


0.99150 


552. 


0.00850 


1.14 ±0.21 


0.830 


2.19 


N6946 


156. 


0.98066 


399. 


0.01934 


1.70 ±0.26 


13.119 


1.46 




126. 


0.98105 


350. 


0.01895 


2.04 ±0.41 


24.494 


1.09 




161. 


0.98857 


461. 


0.01143 


1.43 ±0.24 


13.603 


1.62 


N7331 


136. 


0.99043 


401. 


0.00957 


3 


3.976 






118. 


0.99023 


361. 


0.00977 


3 


6.658 




N7552 


144. 


0.98818 


418. 


0.01182 


1.77 ±0.35 


20.019 


1.43 




115. 


0.99053 


367. 


0.00947 


2.21 ±0.46 


42.249 


1.12 




147. 


0.99277 


483. 


0.00723 


1.47 ±0.26 


21.235 


1.67 


N7793 


176. 


0.96455 


381. 


0.03545 


0.91 ±0.09 


0.869 


2.29 




144. 


0.95400 


334. 


0.04600 


1.02 ±0.20 


1.368 


1.50 




181. 


0.98096 


440. 


0.01904 


0.79 ±0.13 


0.900 


2.62 



'^The parameters /i are the mass fractions of the discrete components at the temperatures Ti. For each 
galaxy, the first Une gives the results of the fits where Ti is constrained in order not to overestimate the 
column density, and the second (third) line where Ti is unrestricted. For galaxies where the S(2) flux is an 
upper limit, only constrained- Ti fits were performed. When OPi?highT < 3, T2 is fixed, and the free-Ti fit 
results are then provided for two different values of T2: 1.14 x r(Sl-S3) (second line) and 1.5 x T(S1-S3) 
(third line). See Section [3T2l for explanations. 

'^The third component is characterized by T3 = 1415 K and /a = 0.00690 (constrained- Ti fit) or T3 = 
1455 K and /a = 0.00477 (free-Ti fit). 

==^3 = 1139K and /a = 0.00088 (constrained-Ti fit), or T3 = 1319K and /a = 0.00019 (frcc-Ti fit). 

^Ta = 1324K and /a = 0.00561 (constrained-Ti fit), or T3 = 1414K and /a = 0.00010 (free-Ti fit). 



Table 5. Masses in the warm and cold phases.^ 



ffala,xv 


M(T > Ti ) 


A^wfcold H9) 
(10^° mol. cm-2) 


M(T > Ti) 
1 MCcoldHol 


CO beam 

1 C4J- V_,0 v_^V_^ ] 


corr. 


map 


ref CO 




N24 


3.13 
















N337 


11.61 


9.4 


0.18 


1521. 


1.52 


IRAC 


E96 




N628 


3.79 


10.5 


0.25 


627. 


1.04 


BIMA 


S05 (B93, Y95) 




N855 


0.50 


2.3 


0.22 


415. 


1.24 


IRAC 


W95 




N925 


1.15 


< 18.7 


> 0.06 


386. 






S05 




N1097 


103.01 


183. 


0.06 


1590. 


1.74 


IRAC 


Y95 (VV98, H93, K03) 




N1266 


12.97 
















N1291 


0.42 


11.7 


0.04 


1452. 


2.33 


IRAC 


T91 




N1316 


0.85 


19.9 


0.006 


1452. 


2.55 


IRAC 


HOI 


1 
1 


N1482 


73.91 


136. 


0.04 


1521. 


1.76 


IRAC 


E96 


Oi 
Oi 

1 


N1512 


6.34 














N1566 


10.22 
















N1705 


0.50 
















N2403 


0.33 


23.8 


0.10 


855. 


0.96 


BIMA 


E96 




N2798 


34.54 


46.5 


0.11 


855. 


2.50 


IRAC 


E96 




N2841 


0.74 


< 27.4 


> 0.03 


370. 






S05 




N2915 


0.0077 


< 6.0 


> 0.02 


1452. 


2.04 


IRAC 


A04 




N2976 


0.28 


10.6 


0.20 


452. 


1.12 


BIMA 


A04 (105, Y95) 




N3031 


0.10 


< 7.9 


> 0.09 


3068. 






S05 




N3049 


3.39 


31.3 


0.02 


573. 


1.48 


IRAC 


C97 




N3184 


2.09 


40.1 


0.07 


441. 


1.23 


BIMA 


S05 (Y95, S93) 




N3190 


11.33 


37.8 


0.09 


380. 


1.14 


IRAC 


L98 




N3198 


3.54 


33.1 


0.10 


415. 


1.31 


IRAC 


B93 




N3265 


6.45 


25.8 


0.06 


227. 


0.76 


IRAC 


G91 




Mrk33^ 


21.35 


16.5 


0.23 


380. 


1.15 


IRAC 


S92 (A04) 




N3351 


11.30 


95.1 


0.05 


606. 


0.85 


BIMA 


P97 (S05) 





Table 5 — Continued 



galaxy 


M{T > Ti) 


Mot (cold H2) 


M{T > Ti) 


CO beam 


corr. 


map 


ref CO 




(106 M0) 


(1020mol.cm-2) 


/ M(coldH2) 


(arcsec^) 








N3521 


4.95 


26.1 


0.22 


201. 


0.68 


BIMA 


NOl (Y95) 


N3621 


1.88 














N3627 


2.39 


196. 


0.01 


271. 


0.95 


BIMA 


S05 (Y95) 


N3773 


2.98 


< 13.5 


> 0.12 


2376. 


5.83 


IRAC 


L05 


N3938 


4.23 


26.4 


0.10 


262. 


0.96 


BIMA 


S05 (Y95) 


N4125 


1.20 


6.5 


0.04 


415. 


1.13 


IRAC 


W95 


N4254 


7.86 


81.8 


0.02 


201. 


0.75 


S99 


NOl (K88) 


N4321 


47.83 


133. 


0.03 


415. 


0.74 


BIMA 


B93 (S05, K96, H93) 


N4450 


5.04 


16.4 


0.07 


131. 


0.90 


BIMA 


S05 


N4536 


102.41 


50.2 


0.11 


855. 


1.08 


S03 


E96 


N4552 


1.91 














N4559'= 


9.29 


13.9 


0.45 


355. 


0.98 


BIMA 


S05 


N4569 


17.48 


200. 


0.02 


118. 


0.73 


BIMA 


S05 (NOl, K88, H93) 


N4579 


2.36 


36.7 


0.01 


118. 


0.42 


BIMA 


S05 (K88) 


N4594 


0.47 


< 6.2 


> 0.04 


346. 






B91 


N4625 


1.85 


9.2 


0.20 


452. 


1.10 


IRAC 


A04 (B03) 


N4631 


8.44 


96.2 


0.10 


1590. 


1.76 


IRAC 


Y95 


N4725 


5.85 


< 54.8 


> 0.03 


210. 






S05 


N4736 


1.20 


81.3 


0.05 


767. 


1.37 


BIMA 


S05 (S98, Y95) 


N4826 


3.00 


250. 


0.04 


1337. 


1.81 


BIMA 


S05 (NOl, Y95) 


N5033 


10.86 


96.0 


0.05 


855. 


1.39 


BIMA 


COl (S05, NOl, E96, A95) 


N5055 


5.98 


122. 


0.06 


645. 


1.23 


BIMA 


S05 (M99, E96, A95, Y95, 


N5194 


1.46 


40.0 


0.04 


474. 


0.98 


BIMA 


S05 


N5194j^ 


0.95 














N5194_B 


3.37 














N5194.C 


6.07 















Table 5 — Continued 



galaxy 


M(r>Ti) iVtot(cold H2) 


M{T>Ti) CO beam corr. map 


ref CO 




(10^ Mq) (102"mol.cni-2) 


/ M(coldH2) (arcsec^) 





iN0iy4_U 


Q C^Q 
0.00 




N5194J] 


3.33 




N5194J 


3.07 




N5194.G 


2.91 




N5195 


5.51 


110. 


N5713 


30.27 


50.8 


N5866 


2.80 


50.9 




0.0047 


7.8 


N6822.B 


0.0062 




N6822.C 


0.0030 




N6946 


11.95 


486. 


N7331 


9.69 


56.0 


N7552 


269.01 


171. 


N7793 


0.10 


25.4 



0.06 


177. 


0.58 


IRAC 


K02 (S89) 


0.08 


1590. 


2.01 


IRAC 


Y95 (Y03) 


0.03 


346. 


1.02 


IRAC 


W03 (T94) 


0.14 


415. 


1.26 


IRAC 


103 


0.03 


816. 


0.95 


BIMA 


S05 


0.07 


346. 


0.94 


BIMA 


199 (NOl) 


0.12 


1452. 


1.86 


IRAC 


A95 (C92) 


0.03 


1452. 


1.87 


IRAC 


195 



^We retain here results from constrained- Ti fits, giving lower masses in the warm phase. Whenever 
the excitation diagrams are ambiguous regarding the value of OPi? high t, we adopt the mass of warm H2 
derived with OPi? high t < 3, which is also smaller than the mass derived with OPi? high t = 3 (see Table HJ. 
References for CO intensities are given in abbreviated form. The quantity "corr." designates the factor 
that was applied to the CO brightness in order to correct for the different apertures of the CO and H2 
observations, and "map" is the image used to derive this factor (see Section |4]). 

'^Mrk33: The H2 mass fraction in the warm phase may be severely overestimated in this blue compact 
galaxy, because of the use of an inappropriate factor to convert CO flux to total H2 mass. Israeli ( 2005 ) 
derive a factor about four times higher than the standard factor used here, which would make the total H2 



mass also four times higher. NGC 6822_A: In their detailed study of this region, Israel et al. ( 2003h derive a 



CO flux to total II2 mass conversion factor twenty times higher than the standard factor, so that the same 
remark as for Mrk33 applies. 

'^NGC 4559 has the highest mass fraction in the warm phase, even when constraining Ti to be near- 
maximal. The warm phase is dominated by gas at ~ 100 K. NGC 4559 is a quiescent galaxy, as obvious from 
the low i^24/-p7.9 flux density ratio within the IRS aperture. The galaxy has little molecular gas, and it is 
possible that it is mostly diffuse, and thus more exposed to FUV radiation. The same remarks may apply, 
to a lesser degree, to other galaxies of the sample such as NGC 628. 
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Table 6. Nuclear X-ray fluxes in the 2-10 keV band and references. 



galaxy 


(10-17 Win-2) 


ref 


note*^ 


N628 


0.59 


K05 


2 


N1097 


173. 


N06 


1 


N1291 


7.45 


K05 


2 


N1316 


10.1 


KF03 


2 


N2841 


1.06 


HFOl 


1 


iNOUOi 


±uzu. 


nr ui 


1 


N3184 


1.12 


K05 


2 


N3627 


< 0.32 


HFOl 


1 


N4125 


1.07 


S04 


1 


N4321 


< 1.25 


HFOl 


1 


N4552 


12.1 


SD05 


1 


N4569 


7.61 


HFOl 


1 


N4579 


264. 


HFOl 


1 


N4594 


120. 


HFOl 


1 


N4725 


7.15 


HFOl 


1 


N4736 


27.4 


SD05 


1 


N4826 


< 1.08 


HFOl 


1 


N5033 


123. 


HFOl 


1 


N5194 


1300. 


S04 


1 


N5195 


4.45 


K05 


2 


N5866 


0.18 


SD05 


1 


N7331 


2.28 


S04 


1 



^1: Quantity given directl y in the c i ted p aper (abbreviated notation as in the list of references). 
2: Quantity computed as in lHo et al.l (|200l[ ) (HFOl) from counts in the 0.2-8 keV band. 




Fig. 1. — Histogram of the equivalent diameters of the areas over which the hne and contin- 
uum fluxes were integrated. The apertures were defined by the intersection of the spectral 
maps in the LH, SH and SL modules of the IRS instrument. The extranuclear regions in 
NGC5194 are shown by the hatched histogram, the regions within dwarfs by the cross- 
hatched histogram, the LINER nuclei by the darker shade and the Sy nuclei by the lighter 
shade. 
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Fig. 2. — Examples of H2 line spectra: the circumnuclear starbursts NGC 1097, NGC 6946 
and NGC 7552, and the three galaxies for which we could estimate the fluxes of higher 
transitions than S(3), NGC 1266, NGC 4569 and NGC 4579. The straight line indicates the 
fitted pseudo-continuum, and the diamond symbols show the wavelength range over which 
the line flux was integrated. 
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Fig. 2. — (continued). 
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Fig. 3. — Critical densities for collisional deexcitation by H2 for the rotational transitions in 
the wavelength range of the IRS instrument. 




Fig. 4. — Examples of how the possible range of apparent O-P-RhighT is determined for 
each galaxy. Thermalization requires the apparent temperatures derived from each pair of 

transitions to be monotonic as a function of upper level energy; in particular, the conditions 
T{S1 — S2) < T{S1 — S3) < T{S2 — 5*3) (indicated by the shaded areas) have to be 
satisfied. The apparent temperatures derived for NGC 3198 arc compatible with OPi? highT = 
3, whereas they require OPi?highT < 3 for NGC 4631 (see text). 
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Fig. 5. — Excitation diagrams. The N^/ ratios are normalized by the S(l) transition. 
The dashed hne indicates the best fit (see text). Whenever OPi?highT < 3 is required by 
the temperature constraints, a dot-dash hne shows for comparison the fit obtained with 
OPi? high T = 3. 
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Fig. 5. — (continued). 
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Fig. 5. — (continued). 
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Fig. 5. — (continued). 
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Fig. 6. — Left: Histograms of the derived column densities in fits where the lower tem- 
perature is constrained (top) and where it is free to vary (bottom). In the latter case, the 
solid line indicates results with low T2 values and the dashed line with high T2 values, in 
OP-RhighT < 3 fits (see text). Righ t: Column densi ties as a function of the total stellar mass 
of the host galaxy, estimated as by lLee et al.l (l2006l ). Star- forming nuclei are represented as 
diamonds, the extranuclear regions in NGC 5194 and NGC 6822 as squares, LINER nuclei as 
star symbols and Sy nuclei as crosses. By definition, dwarf galaxies have M^^ost < 10^ '' M^r 
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Fig. 7. — Fraction of H2 in the warm phase (T > Ti) as a function of Ti, the lower tem- 
perature of the two components fitted to the rotational lines. Nuclei classified as purely 
star-forming are shown as diamonds, and regions in dwarf galaxies as filled diamonds; nuclei 
classified as LINER or Sy are shown as star symbols or triangles, respectively. The symbol 
size is increased for targets with the most robust CO flux estimates (see text) . In the free-Ti 
fits with OPR highT < 3, results with the two values of T2 considered here (see text) are con- 
nected by a line segment; the fits with T2 — 1.5 x ^(Sl — S3) produce higher temperatures 
and lower warm H2 masses than the fits with T2 = 1.14 x T{S1 — S3). 
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Fig. 8. — Ratio of the power emitted in the sum of the S(0) to S(2) transitions to: (a) the 
total infrared power; (b) the [Sill] hne power; (c) the power emitted in the aromatic bands 
within the IRAC4 filter (assuming a filter width of 13.9 THz and subtracting the stellar 
emission as mentioned in Sect. 12.21) : (d) the power emitted at 24 /im within the MIPSl 
filter (assuming a filter width of 3.1 THz). The symbol coding is as in Fig. [71 Excluding 
the regions within NGC6822 (the targets with the smallest projected aperture), as well as 
NGC1705 and NGC2915 (the galaxies with the smallest H2 brightness), the average and 
dispersion of each logarithmic power ratio is: (a) —3.36 ± 0.15 for HII nuclei and complexes; 
and -3.12 ± 0.24 for LINER and Sy nuclei; (b) -0.53 ± 0.17 and -0.04 ± 0.42, respectively; 
(c) -2.19 ±0.10 and -1.80 ±0.34; (d) -1.85 ± 0.28 and -1.51 ± 0.31. Note that the H2 flux 
fabscissa) is almost nronortional to the Hq brightness, because of the nuasi-uniform beam. 
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Fig. 9. — Power ratios of H2 to aromatic bands (shown in Fig. [8t) and 24 /im emission to 
aromatic bands as a function of P24 = 1.05((z/24F24 — 0.14z/7 QF7 C)Hiist)/ ('^7i-p7i + ^^i.'^g - F'i.^fi) — 
0.035)° ''^, a quantity closely related to /(Ldust ; t^ad > 100) defined by lDraine et al.l (120071 ) . 
which is the fraction of the total dust luminosity coming from regions with radiation field 
intensities more than 100 times the average local field. Only HII nuclei and complexes 
are shown. Dwarfs are represented by filled diamonds and the extranuclear regions within 
NGC5194 by thick squares. 
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Fig. 10. — Ortho to para ratios as a function of the surface brightness in the sum of the 
S(0) to S(2) transitions and as a function of the surface brightness in the 24 /im band. For 
a number of sources, based on the excitation diagram, we adopted a fixed OPi? high t = 3 
(see text). Galaxies with upper limits in the S(2) line, or with indications of possibly non- 
neghgible optical depth at 10 /im, are shown as lower limits of OPi?highT- For free-Ti fits, 
results obtained from the two values of T2 considered here are connected by a line segment. 
The symbol coding is as in Fig. U\ except that LINER and Sy nuclei are both shown as star 
symbols. 
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Fig. 11. — Excess H2 luminosity as a function of 2-10 keV X-ray luminosity. The excess 
H2 emission is defined from the relation between H2 and aromatic band power, shown in 
Figure [Ht, as the difference between the total H2 emission and the quantity 10~^'^^ x -Fy gdust, 
which defines the upper envelope of HII nuclei. Galaxies with no H2 excess according to this 
definition were arbitrarily placed at an ordinate of 2.25 below the dashed line. The size of 
the symbols is proportional to the fraction of excess H2 to total H2 emission, coded as in 
Fig. [T2I Sy nuclei are marked by an overlying horizontal bar, and star-forming nuclei by a 
square. The large scatter and the high ratios of excess H2 to X-ray luminosities argue against 
X-rays playing a dominant role in the additional (nonstellar) excitation of H2 in LINER/Sy 
nuclei (see text). 
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Fig. 12. — Heating efficiency required to account for the excess H2 emission by supernova 
remnant shocks (ratio of the power in the sum of the S(0)-S(2) hues to the total mechanical 
power), as a function of the excess H2 luminosity, defined as in Fig. [HI The dashed line 
indicates the maximum efficiency expected if all the mechanical power pr oduced by supernova 
remna nts is absorbed in molecular clouds, for the model parameters of iKaufman fc Neufeld 
( 119961 ). The size of the symbols is proportional to the fraction of excess H2 to total H2 
emission. Sy nuclei are marked by an overlying horizontal bar. The regions B and C within 
NGC 6822 are the only non-AGN targets in this figure. 



